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Abstract 

t-h We review the evolution of the studies of diffractive processes in 

the strong interaction over the last 60 years. First, we briefly outline 
the early developments of the theory based on analyticity and unitar- 

fT) ity of the S'-matrix, including the derivation and exploration of the 

Regge trajectories and related moving cuts. Special attention is paid 
to the concept of the Pomeron trajectory introduced for description 

Z^pr of total, elastic and diffractive cross sections at high energies and to 

T — I the emergence of the dynamics of multi-Pomeron interactions. The 

^ role of large longitudinal distances and color coherent phenomena for 

the understanding of inelastic diffraction in hadron-hadron scattering 
and deep inelastic scattering is emphasized. The connection of these 
C3 phenomena to the cancellation of the contribution of the Glauber ap- 

proximation in hadron-nucleus collisions and to the understanding 
of the Gribov-Glauber approximation is explained. The presence of 
different scales in perturbative QCD due to masses of heavy quarks 
has led to the emergence of numerous new phenomena including non- 
universality of the slopes of Regge trajectories made of light and heavy 
quarks and non-universal energy dependence of elastic cross sections. 
The application of the perturbative QCD techniques allowed us to 
calculate from the first principles the interaction of small transverse 



size color singlets with hadrons leading to the development of the 
quantitative theory of hard exclusive reactions and to the successful 
prediction of many regularities in hard large mass diffraction. It also 
led to the prediction of the phenomenon of complete transparency of 
nuclear matter in QCD in special processes. The conflict of pertur- 
bative QCD with probability conservation for high energy processes 
of virtual photon-nucleon scattering is explained. Some properties of 
the new QCD regime are outlined. 

1 Introduction 

The aim of this chapter is to demonstrate that the phenomenon of diffrac- 
tion in quantum chromodynamics (QCD) is a formative playground for the 
fundamental ideas and methods of theoretical physics. The deep disappoint- 
ment of scientific society in the quantum field theory paradigm formulated 
by L. D. Landau at the Rochester Kiev conference (1959) pQ was based on 
the zero-charge problem in pre-QCD quantum field theories. As a result, 
the concept of the scattering matrix (S-matrix), where all quantities are in 
principle observable, became popular and replaced studies within the quan- 
tum field theory framework. The idea was that unitarity of the S-matrix, 
its analytic properties and exact symmetries will allow one to avoid dealing 
with point-like interactions characteristic for a quantum field theory such 
that the need for the ugly procedure of renormalization with all its puzzles 
will disappear. This approach led to the development of such new concepts 
as single and double dispersion relations, Regge trajectories, the Pomeron 
calculus, string models, etc and to the prediction of new phenomena. The 
discovery of asymptotic freedom in QCD in the late sixties to early seven- 
ties justified the space-time description of high energy processes in QCD, 
which is absent within the S'-matrix concept. The account of the space-time 
evolution of high energy processes allowed one to predict a variety of strik- 
ing new QCD phenomena such as color fluctuations, complete transparency 
of nuclear matter under special kinematic conditions, formation of a new 
QCD regime of the maximally strong interaction, etc., all of which are ab- 
sent in the S'-matrix theory. In the first part of this chapter, we consider 
the phenomenon of diffraction in the S-matrix theory and then discuss new 
diffractive phenomena that emerge in QCD studies. 

The basic ideas of the S-matrix approach are unitarity of the S-matrix 



in all physical channels and analyticity of scattering amplitudes in the com- 
plex planes of energies and momentum transfers that leads to analyticity of 
amplitudes in the plane of the angular momentum. The S'-matrix approach 
justifies the concept of Regge trajectories. The assumption of the domi- 
nance of amplitudes of high energy processes by the Pomeron trajectory 
exchange predicts an increase with energy in the radius of the interacting 
hadron (shrinking with energy of the forward peak in two-body exclusive 
processes) and, therefore, the dominance of peripheral collisions, the univer- 
sal dependence on energy of the total and elastic cross sections and cross 
sections of diffractive processes. The very existence of Pomeron moving cuts 
follows from unitarity of the S'-matrix in the crossed channel. The prediction 
and experimental discovery of the large mass (M 2 ^> m 2 ) triple Pomeron 
diffraction proves the non-zero value of the effective triple Pomeron interac- 
tion. 

The modeling of the contribution of Pomeron moving cuts found black- 
ening of interactions at central impact parameters since the contribution of 
the single Pomeron exchange grows with energy. However, the fraction of 
the total cross section due to the elastic scattering slowly grows with energy. 
For the current LHC energy: 

Elastic (PP)/<?tot(PP) « 0.25 . (1) 

Thus, the pp interaction is still very far from the regime of complete ab- 
sorption where this ratio should be close to 0.5. Note, however, that for the 
central pp collisions, almost complete absorption has been observed at Fermi 
National Accelerator Laboratory (FNAL) and the LHC, which corresponds 
to partial amplitudes being close to unity. Thus an energetic proton, when 
interacting with the proton target, behaves as a grey disc with a black spot 
in the center. The size of the black spot rapidly increases with energy. 

The assumption that the amplitudes of high energy hadron-hadron colli- 
sions depend on one scale was challenged by the discovery of J/ip, T mesons 
- bound states of heavy quarks: cc and bb. The radii of these quarkonia 
states are significantly smaller than for hadrons made of light quarks. The 
interaction of QQ quarkonia with hadrons made of light quarks is decreasing 
with rriQ in the non-perturbative and perturbative QCD domains. 

Analyses of the ladder diagrams for cross sections of deep inelastic scat- 
tering off a hadron target found that longitudinal distances dominating in 
the scattering process are linearly increasing with energy in QCD. As a re- 
sult, at sufficiently large energies they exceed by far the length of the target, 



the transitions between different configurations in the projectile slow down 
and the interaction can be described as a superposition of the interaction of 
instant quark-gluon configurations within the projectile. This feature leads 
to the fluctuations of strengths of hadron-hadron and hadron-nucleus inter- 
actions and to the exact cancellation of the Glauber model contribution to 
hadron (nucleus)-nucleus collisions. This cancellation follows directly from 
the analytic properties of amplitudes and/or energy-momentum conserva- 
tion. The Gribov-Glauber model replaces the Glauber model in high energy 
processes where diffraction is a shadow of inelastic processes. 

The fluctuations of strengths of the interaction in hadron-hadron col- 
lisions found an explanation in QCD as being due to the color screening 
phenomenon. The fluctuations of strengths of the interaction within the vir- 
tual photon wave function have been observed directly in the significant cross 
section of leading twist diffraction in deep inelastic scattering (DIS), which 
is predicted to be negligible in perturbative QCD (pQCD) because of the 
absence of free quarks and gluons. However, it was observed in ep DIS at 
HERA that the cross section of diffractive processes constitutes ~ 10% of the 
total cross section at x ss 10~ 3 and the ratio of the diffraction cross section 
to the total one is practically energy, and Q 2 - independent (Q 2 is the photon 
virtuality) . 

QCD dynamics predicts the existence of hard diffractive phenomena that 
are higher twist effects. The significant difference in momentum scales char- 
acterizing hard and soft (non-perturbative) processes allows one to prove the 
factorization of hard processes from soft ones and to calculate cross sections 
of hard diffractive processes. The processes of elastic photoproduction of 
mesons with hidden heavy flavors off the proton target and elastic electro- 
production of light mesons observed at HERA revealed an interplay of the 
dependence of the cross sections on energy, the photon virtuality Q 2 , and the 
momentum transfer t that is close to that expected in pQCD. 

Diffraction in high energy processes is a shadow of inelastic processes so 
that, for sufficiently small x, an increase with energy of the structure func- 
tions of nucleons and nuclei, which is predicted in pQCD approximations, 
runs into conflict with the probability conservation at small impact parame- 
ters. The range of central impact parameters, where the regime of complete 
absorption dominates, increases with an increase of energy (for fixed Q 2 ). In 
this regime at ultrahigh energies, a(pp) oc ln 2 (s/s ) and a{^*p) oc ln 3 (s/s ). 
The complete absorption regime is possibly reached for the gluon distribu- 
tion at the central impact parameters and at Q 2 of the order of a few GeV 2 



in electron-proton collisions at HERA. Physics related to the formation of 
strong gluon fields at sufficiently small x ~ Q 2 /s and small impact parame- 
ters can be probed at the LHC and Large Hadron-Electron Collider (LHeC). 
In this new QCD regime the expansion over powers of 1/Q 2 (twists) becomes 
meaningless. Also, there arises the question whether continuos symmetries 
such as conformal and scale invariances characterizing pQCD and new QCD 
regimes are different. 

No significant violation of the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi 
(DGLAP) approximation for the structure functions integrated over impact 
parameters is predicted at achieved x. This is because the structure func- 
tions are dominated by the scattering at the impact parameters growing 
with energy where the interaction remains weak even though it reaches the 
black limit at the central impact parameters. Moreover, no noticeable slow- 
ing down of an increase of the structure functions with energy is expected 
at even smaller x since the black limit contribution rapidly increases with a 
decrease in x at the energies achievable in laboratory. 

It was suggested to sum the leading a s ln(xo/x) terms in the kinematics of 
fixed Q 2 and x — > - leading logarithmic (LL) approximation. There were 
derived formulas for the collision of two small size (~ 1/Q) wave packets, 
j*(Q 2 ) + 7*(Q 2 ) — > X, in which case the diffusion in the k t space is sup- 
pressed (which may work for limited range of energies) . The same formulas 
are often applied to the scattering of a small size (~ 1/Q) wave packet (7*) 
off the proton target that has a size of ~ 1/(2???^). In this case the neglected 
within the LL approximation diffusion to small parton momenta within the 
parton ladder is rather important. Energy and momentum of the final states 
calculated within the leading log approaches are significantly different from 
that for initial state especially within the leading a s ln(xo/a;) approximation. 
This violation follows from the choice of the kinematical domain characteris- 
tic for LL approximation. Conservation of energy-momentum is guaranteed 
after resummation over series of LO, NLO, NNLO... approximations. This 
property of approximation explains negative sign and huge value of next-to- 
leading order (NLO) "corrections" and requires the development of resum- 
mation approaches. However, diffractive processes were not considered yet 
in the resummation approaches. 

The chapter is organized as follows. In section [2] we define the kinemat- 
ics characteristic for diffractive phenomena in hadronic collisions and briefly 
review the ^-matrix approach: analytic properties of amplitudes of high en- 
ergy processes in the energy, momentum transfer and orbital momentum 



planes. We also explain how the concept of Regge trajectories in the angular 
momentum plane arises in the relativistic theory of the scattering matrix. 

In section [3j we remind basic properties of Regge trajectories and explain 
that the linearity and universality of Regge pole trajectories is confirmed by 
comparison of Regge trajectories with the data on hadron resonances made 
of light u, d, and s quarks. This linearity and the assumption that the 
amplitude is dominated by the Regge pole trajectory exchange allows one to 
reproduce the observed dependence of cross sections of two-body processes 
with non-vacuum quantum numbers in the crossed channel on energy. We 
also explain that the assumption that all Regge trajectories have the same 
universal slope does not hold for the trajectories with hidden and open heavy 
flavors. It contradicts the quarkonium models and the data. We also explain 
how moving cuts accompanying the Regge pole follow from the unitarity of 
the S'-matrix in the plane of angular momentum. 

In section |4| we consider the hypothesis that the dominant exchange in 
the amplitudes of high energy processes, is the one by the Pomeron trajectory 
accompanied by Pomeron moving cuts — Pomeron calculus. The dominance 
of peripheral collisions and the important role of Gribov diffusion in the im- 
pact parameter space are explained. A brief comparison with data shows that 
such predicted basic features of high energy processes as the universal energy 
dependence of all high energy processes and the shrinking with energy of the 
diffractive peak agree with the data. We discuss properties of multi-Pomeron 
interactions, the evidence for triple Pomeron interactions and their role in 
diffraction at the Tevatron and LHC The impact of multiple rescatterings 
of Pomerons on the elastic differential cross section is also briefly discussed. 

In the framework of the Pomeron calculus, we also explain an onset of 
the regime of the complete absorption at small impact parameters, some of 
its properties and compare it briefly with the selected FNAL and LHC data. 
Implications for the value of the slope of the Pomeron trajectory in the regime 
of complete absorption at small impact parameters are briefly discussed. 

In section [5] the space-time evolution of high energy processes and the 
linear increase with energy in longitudinal distances in the scattering process 
are discussed. As a result, the contribution of planar diagrams — known as 
the Glauber model for hadron-nucleus collisions — is cancelled out. Moreover 
the contribution of the planar diagrams (the Glauber approximation) violates 
energy-momentum conservation. 

At the same time, the contribution of non-planar diagrams can be rewrit- 
ten in the form of the Glauber approximation but with an additional inelastic 
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shadowing term — the Gribov-Glauber model. 

In section[6]we explain that in QCD, the increase of longitudinal distances 
with the collision energy leads to a variety of coherent phenomena which we 
refer to as the color fluctuation phenomena. They include the presence of 
inelastic diffraction at the zero angle and the processes where hadrons fluc- 
tuate into small - size configurations and interact with the small strength so 
that nuclei do not absorb them — color transparency (CT). CT allows one to 
prove QCD factorization theorems for a number of processes. At the same 
time, the concept of color fluctuations allows us to bridge the gap between 
the fluctuations to small- and large-size configurations and to reconstruct the 
distribution over cross sections for projectile hadrons and photons. It also 
gives us a physically transparent interpretation of the Gribov-Glauber model 
that could be applied to modeling proton (nucleus)-nucleus collisions. The 
concept of color fluctuations allows us to build the QCD-improved aligned 
jet model. The significant cross section of diffraction in deep inelastic small 
x processes observed at HERA, its Q 2 and energy dependencies are direct 
confirmation of the important role of color fluctuations in high energy pro- 
cesses. We also briefly review the concept of the perturbative Pomeron in 
pQCD. 

In section [7| we explain the QCD factorization theorem for hard exclusive 
processes, derive basic characteristics of hard diffractive processes and briefly 
compare the derived formulas with the data obtained at FNAL and HERA. 
We point out that complete transparency of nuclear matter for special hard 
diffractive processes has been predicted and confirmed by the FNAL data. 

In section |8j we discuss the onset of the new QCD regime and its basic 
features in the limit of fixed Q 2 and x — > 0. The conflict between pQCD 
calculations and probability conservation in the collisions at central impact 
parameters, the onset of the black disc regime (BDR), and competition be- 
tween the soft QCD and pQCD contributions are explained. 

Conclusions are presented in section [9] . 



2 General properties of the scattering ampli- 
tude 

2.1 Kinematics 

We consider first the scattering of two particles: 

a + b^ a + b. (2) 

The amplitude of this process depends on the four-momenta of the col- 
liding particles pj. An account of energy-momentum conservation gives the 
constraint: 

Pa + Pb = p' a + Pb ■ ( 3 ) 

Lorentz invariance restricts the number of independent variables. The 
convenient variables are the square of the energy in the s-channel center of 
mass: s = (p a + Pb) 2 , and the squares of the momentum transfer between a 
and a' t = (p a — p' a ) 2 , and between a and b': u = (p a — p' b ) 2 . These three 
variables are not independent. Energy-momentum conservation (Eq.|3]) leads 
to 

s + t + u = 2[m 2 a + m 2 b ] . (4) 

If the scattered particles are the lowest-mass states in the channels with given 
quantum numbers, there exist three physical channels where the scattering 
process is allowed: s > (m a +m ) 2 , (t,u) < 0; t > max{4m^,4m^}, (s,u) < 0; 
u > (m a + rrib) 2 , (t,s) < 0. Thus the physical meaning of the variables s, 
t and u is that each of them is equal to square of the center-of-mass energy 
of colliding particles in the corresponding center of mass of the physical 
channel: a and b in the s-channel, a and a in the t-channel and a and b in the 
w-channel. The amplitudes in all three channels are interrelated by rotation 
of the four-momenta of the particles. 

2.2 S-matrix approach 

Before the advent of quarks and later QCD the most important ideas of the 
theory of strong interactions were suggested within the concept of the scat- 
tering matrix, 5*. The matrix elements of the S'-matrix describe amplitudes 
of the scattering processes where hadrons in the initial and final states are 



outside the interaction region. In the physical region of any physical process, 
the S- matrix is restricted by its unitarity, i.e., by probability conservation: 

SS j = 1. (5) 

To single out the contribution when no interaction occurs, the T-matrix is 
introduced: S = 1 + IT. The unitarity condition for the T-matrix has the 
following form in the s-channel: 



Im T(a + b ->• a + b) = - / ^ (a + b\ T \ri) dr n (n\ T\a + b) ] , 



(6) 



where dr n is the phase volume for the state n. The above equation can be 
rewritten as the optical theorem which relates T matrix with the total cross 
section: 

ImT(a + b ->■ a + b) = sa tot (a + b) , (7) 

where s = (p a +Pb) 2 - 

The hope was that in a relativistic theory, the conservation of probability, 
i.e., the S'-matrix unitarity (Eq. [5]), and threshold singularities in the crossed 
channel would substitute the non-relativistic concept of the potential. 

2.3 Brief summary of analytic properties of amplitudes 
in energy, momentum and orbital momentum planes 

The aim of this subsection is to briefly remind of basic ideas and results 
obtained in the S'-matrix approach. 

One starts with imposing causality in the form of the Lehman, Symanzic 
and Zimmermann (LSZ) representation of the amplitude as the Fourier trans- 
form of the matrix elements of the retarded commutator of currents [2] . This 
leads to the assumption that the amplitudes of the physical processes are the 
boundary values of the same analytic function of the energy and momentum 
transfers. The singularities of the amplitudes are given by the thresholds for 
the physical processes [3], see the discussion below. 

S. Mandelstam proposed the double dispersion representation which takes 
into account the singularities both in the momentum transfer and in the 
energy planes [1]: 

4,, 'f^_4f) + 'f^ "«<<«'><'> - 



71-2 Am* (s'-s){t'-t) TT 2 J Am l (u'-u)(t'-t) 
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71-2 Am 2 ( S ' — S )( U ' ~ U ) ' 

where we take m a = mb for simplicity. All denominators are understood as 
having imaginary parts: (V — s) -1 = (s' — s — ze) _1 . (This condition selects 
the outgoing wave.) In the following we will not need the exact form of the 
spectral densities pij which were supposed to follow from unitarity of the 
S- matrix and its analytic properties [5]. 

Dispersion representations over the variables s and t follow directly from 
the double dispersion representation. Also, they can be derived directly from 
the theoretical analysis of the LSZ representation of the scattering amplitude, 
see the discussion and references in [0] . The dispersion representation of the 
scattering amplitude in the energy plane is: 

A/ , 1 f°° , Jm B A(8 t ,t) 1 r° , ,Im u A(u',t) 

A(s,t) = - / ds v ' ; + - / dv! v ' ; + subtractions . (9) 

^ h™? s ' - s n J 4m 2 v! -u 

By definition subtractions do not have imaginary part in variables s and u. 
It follows from unitarity of the S'-matrix and analyticity that in the physi- 
cal region, the scattering amplitudes at large energies are restricted by the 
condition: ImA < csln (s/s ) |7J. For the amplitude symmetric under the 
transformation s — > u, the subtraction term is constant. For the amplitude 
antisymmetric under the transformation s — Y u (negative signature), the 
subtraction term is oc s. 

The dispersion representation over the momentum transfer t has a similar 
form: 

A , . 1 r° , ,Im u A(u',t) 1 I" 00 Jf Im t A(t',s) 
A(s, t) = — / du 1 — / dt h subtractions . 

7T J im ? W —U 71 J 4m 2 V -t 

(10) 
Emergence of the concept of Regge poles in the relativistic theory was enabled 
by the combination of Eq. [8] for the analytic continuation of the decomposi- 
tion of the scattering amplitudes over partial waves to the crossed channel 
with the S'-matrix unitarity condition in the crossed t-channel. 

One starts with the observation that the S'-matrix unitarity condition 
becomes diagonal if the conservation of the angular momentum is taken into 
account. Hence it is convenient to decompose the amplitude over partial 
waves: _ 

A{s, t) = 87iJ2 fi( s W + l)Pi(z) , (11) 
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where z = 1 + 2t/(s — 4m 2 ); Pi(z) are the Legendre polynomials. For sim- 
plicity, we consider collisions of hadrons with spin zero. (A generalization to 
the case of scattering of particles with non zero spin is straightforward but 
would make formulae unnecessary lengthy.) 

m = I J\ zPl ( Z )^ . (12) 

The normalization of fi is chosen so that unitarity of the S'-matrix has the 
form: 

hamfi(s) = ~fi{s)f*{s) + positive terms. (13) 

In high energy processes, orbital momenta essential in the scattering pro- 
cess are large. Hence it is legitimate to substitute the sum over the orbital 
momenta I by the integral over the impact parameters b, 1 + 1/2 — pb, where 
p is the center of mass momentum. (The factor of 1/2 follows from the 
necessity to reproduce the formulae of the semi-classical approximation in 
non-relativistic quantum mechanics.) At large energies, p ~ \fs/2. Thus we 
derive the impact parameter representation of the amplitude: 

A(s,t) =4tts fdbbf(b,s)J (qtb) = (2s) f d 2 bexp(tq t •&)/(&, s) , (14) 

where q t = psin(9) and 9 is the cm. scattering angle. In the derivation 
we used the asymptotic expression for the Legendre polynomials at large /: 
Pi{9) rs Jq((1 + 1/2)9) and the integral representation of the Bessel function 
Jo- Jo(q) = (l/27r)/ o 7r exp(igcos(0))rf0. 

We can use the dispersion representation of the amplitude over the mo- 
mentum transfer and properties of the Legendre functions of the second kind 
to derive the representation that can be easily continued into the complex 
plane of the orbital momentum: 

i / QM^l dz . (15) 



2ni J c 8tt 

Here the contour of integration encircles the [—1, 1] interval on the real axis. 
The integrand has singularities outside the contour at 

Zl = l-((4m 2 )/(s-4m 2 )), z 2 = -1 + ((4m 2 )/(s - 4m 2 )) , (16) 

corresponding to the singularities of the amplitude at t — 4m 2 and u = 4m 2 . 
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The Legendre functions of the second kind satisfy the relation: 



1 '! 



Qi(z) = -I P l {z')dz'l{z'-z). (17) 



The advantage of this function is that for z 3> 1: 



Qi(*)«i, (18) 

and that for — 1 < z < 1, 

Q,(« + ie)-Q,(«-ie) = -»7rfi(z). (19) 

For the discussion in the next section, it is important to derive the rep- 
resentation of the partial wave in t-channel. For sufficiently large /, the 
integration contour in Eq. [T5] can be deformed around the singularities of 
the amplitude since the integral over the large circle is equal to zero. Thus, 
another representation arises: 

^ar ^'^^ Hfi f «'^»' t . (20) 

n J Z1 8vr ix J Z2 8vr 

The presence of the factor of (—1)' = exp(iirl) precludes analytic con- 
tinuation of the amplitude to the complex plane of the orbital momentum 
since the factor (—1)' = exp(iirl) increases rapidly with Iml. To remove the 
factor of (— 1)', it is convenient to introduce the functions f t (t) which are 
symmetric and antisymmetric under the s — > u transposition, respectively — 
functions with the positive and negative signature. Thus, 

rim- 1 rQi(z)A t (z,s) ^ 1 [°° Qi{z)A t {z,u) 



ft(t) = - ^ ' " - J ±- "" 'T , (21) 



7T J Z1 87T 7T J Z2 <Sn 



which provides the analytic continuation of the partial waves to the complex 
plane of the angular momentum. Above formulae are known as the Gribov- 
Froissart projection [8]. 

3 Regge poles in the /S-matrix theory 

3.1 Regge poles and ^-channel unitarity 

T. Regge found in the non-relativistic quantum mechanics that the scat- 
tering amplitude in the unphysical region corresponding to large imaginary 

12 



scattering angles, cos(#) — > oo, has the following form: 

A(s,t) (x cos(6) l{E) , (22) 

where 1(E) is the eigenvalue of the operator of the orbital momentum I at a 
given energy E. The eigenvalues of the energy E n follow from the condition 
that / = n, where n is an integer number. Thus the concept of the Regge 
trajectory 1(E) allows one to describe both the energy eigenstates and the 
asymptotic behavior of the amplitude [9]. In a relativistic theory, the s- 
channel cos(#) oc s and E should be substituted by t. 

S. Mandelstam observed that in a relativistic theory, the kinematics of 
large s and fixed t<s corresponds to physical processes with s ^> 4m 2 which 
are usually called the crossed channel with respect to the t-channel processes. 
As cited in [TO], he suggested that the Regge pole behavior, would allow for 
a simple description of bound states. The key tool for the derivation of the 
basic properties of Regge trajectories and the calculation of the amplitudes of 
high energy processes is the partial amplitudes in the t-channel analytically 
continued to the complex plane of orbital momentum — the Gribov-Froissart 
projection discussed in the previous section. Regge trajectories describe the 
sum of poles of these amplitudes which follow from unitarity of the ^-matrix 
in the crossed channel. The concept of the Regge trajectory l(t) is useful for 
the description of hadron resonances with the same quantum numbers (except 
for spin) and for the calculation of amplitudes of high energy processes. 

Important properties of Regge trajectories follow from the two-particle 
t-channel unitarity condition continued to the angular momentum plane: 

\(f?(t + ie)) - ft(t - ie)) = (l/2)/±(* + ie)ft(t - ie) , (23) 

which can be rewritten as 

^(*) "/?=*(*) = /( ± (*)/ ± *(*)- ( 24 ) 

In the proof one uses the observation that for t < 4m 2 , / are real and hence 
f ± (t — ie) = [/ ± (t + ie)]*. The two-particle t-channel unitarity condition is 
exact for 4m 2 < t < 16m 2 and allows one to prove the existence of Regge 
poles in relativistic amplitudes and to establish some of their properties [Til 

The Gribov-Froissart projection has pole at I = l(t): 

f l (t) = c/(l-l(t)). (25) 
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Taking into account the real and imaginary parts of the trajectory /(£), it 
is easy to find out that in the vicinity of I, the amplitude is described by 
the Breit-Wigner formulas, see also [13J. Thus the concept of the Regge 
trajectory (which is often called the moving trajectory) describes hadronic 
resonances with the same quantum numbers, except for the spin. 

Using these equations and iterating one Regge trajectory, it is easy to 
show that the Regge trajectory generates the moving pole singularities in 
the complex plane of the angular momentum / [121 H3] • 

3.2 Regge poles and high energy behavior of ampli- 
tudes of physical processes 

To demonstrate the role of Regge poles in high energy processes in the rela- 
tivists theory, it is convenient to use the method applied by Sommerfeld to 
the problem of diffraction of radio waves around Earth. The task is to find 
an analytic function of I which coincides with fi(t) for integer points 1 = 0, 
1,2,.... 

The decomposition of the positive signature amplitude over partial waves 
in the t-channel diverges with an increase in s because the Legendre polynomials- 
being the functions of z = 1 + 2s/ (t — 4m 2 ) — increase as powers of s: 

exp (m/2)Q 
P i (cosh(a)) / ^ 00 oc === . (26) 

y 27rsin(a) 

The partial amplitudes are restricted by the S-matrix unitarity and their 
imaginary parts are positive. Thus one needs to continue partial waves an- 
alytically to the angular momentum plane. The procedure was explained in 
the previous section. 

The first step is to identically represent the amplitude as a contour inte- 
gral over the orbital momenta I > around the real axis: 

A+(s, t) = (1/i) / dltf(t)f+(t)(2l + 1)P,(1 + 2s/(t - Am 2 )) , (27) 

Jc 

where £ z + (£) is called the positive signature factor. It can be written as 

£+(£) = [l + (-1)**']/ sin(TrZ) = exp z(tt//2)/ sin(7r//2) . (28) 



Taking the residues over the poles of l/sin(7r//2) would recover Eq. 11 
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For the amplitude antisymmetric with respect to the transposition s — > u, 
the signature factor is: 

f_(0 = expi(7r//2)/cos(7r//2) . (29) 

It follows from the location of singularities in the plane of t that fi(t) oc 
exp(— jjlI) for I — > oo, where /i is the minimal mass in the singularities over t. 
Hence, for the integration contour C, one can take the straight line between 
the points Iq — zoo and Iq + ioo. The contour can be moved to the left until 
it encounters singularities of the amplitudes in I that do not allow further 
shifting of the counter to the left. In this discussion, l Q > —2 since the positive 
signature amplitude cannot decrease with s faster than 1/s 2 . This property 
follows from the fact that imaginary part of this amplitude is always positive. 
In addition, one should take into account poles and cuts in the /-plane. Thus, 
the expression for the positive signature amplitude reads: 

t-lo+oo 

A + (s,t) = (l/i) / die(l)lt(t)(2l + 1)P/(1 + 2s/(t - 4m 2 )) + A , (30) 

where A is the contribution of the Regge poles and moving cuts. A similar 
analysis can be performed for the negative signature amplitude. 

The contribution of the Regge pole fi(t) oc 1/(1 — a(t)) to the amplitude 
A has the form: 

9a,a'(t)g bM (t)(s/s r^at)- (31) 

Factorization of the dependence of the amplitude on the properties of the 
particles a and b follows from the unitarity condition for the partial waves in 
the t-channel [T2] , 

V. Gribov has demonstrated [TEl] that the textbook models that assume 
that the total cross section of hadronic collisions at large s is energy inde- 
pendent and the t dependence of the elastic amplitude does not depend on s, 
i.e., A(s,t) = isf(t), are incompatible with probability conservation in the 



crossed channel as given by Eq. 24 Such a behavior corresponds to the fixed 
pole in the orbital momentum plane. However, a single pole 1/(1 — 1) on the 
left-hand side of this equation cannot be equal to a double pole 1/(1 — l) 2 on 
the right-hand side of it. 

Common wisdom based on the S'-matrix approach was that any hadron 
is a bound state of other hadrons but not of some elementary constituents. 
This concept was implemented by assuming that hadrons belong to Regge 
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trajectories. In quantum chromodynamics, hadrons are bound states of el- 
ementary particles — quarks and gluons — and physical states contain no free 
quarks and gluons (the hypothesis of confinement of quarks and gluons). 
This makes the Regge trajectory description of hadrons even more plausible. 

3.3 Regge trajectories 

G. Chew and S. Frautchi suggested to describe the spectrum of hadrons in 
terms of the Regge pole trajectories |13j . The mass of a hadron follows from 
the equation: 

a(t h ) = J h , (32) 

where Jh is the spin of the resonance and th its mass. The scattering ampli- 
tude has poles for these values of I and t. 

The observed spectrum of hadrons made of light quarks as well as cross 
sections of exclusive processes with non-vacuum quantum numbers in the 
crossed channel are well described by the Regge trajectories linear in the 
momentum transfer t [13] : 

a(t) = a + oc't . (33) 

The data prefer practically the same slope a' for all hadrons made of the 
light u, d, and s quarks: 

a' « 1 GeV~ 2 . (34) 

The values of the intercept ao for the leading trajectories are the following: 
a p (t = 0)«0.5, O4 2 (£ = 0) ^0.3, o*(t = 0)«0. (35) 

For completeness, we also enumerate other important results which, how- 
ever, we will not use in this chapter: 

(i) An analysis of the experimental data on meson resonances and exclu- 
sive processes indicates that the meson trajectories with positive and negative 
signatures are close, see [H] for a review and references. 

(ii) Fermion trajectories should also contain the term ex y/t, otherwise 
the trajectories with opposite signatures will be degenerate [IT]. However, 
the data on two body processes dominated by fermion exchanges are very 
limited. 

(iii) There has been proposed a technique of calculation of the correspond- 
ing QED amplitudes in the angular momentum plane [18] which allows to 

16 



establish whether constituents of QED belong to moving Regge pole trajec- 
tories. A more powerful technique of identifying and calculating the leading 
Feynman diagrams in a quantum field theory containing vector particles has 
been developed in [T7] . 

(iv) It has been shown that in non-abelian gauge theories massive vector 
mesons [19] and fermions [20J are reggeized in the perturbative regime. In 
QCD the presence of infrared singularities related to the zero mass of the 
gluon requires to take into account addition infrared factor. 

(v) It has been observed that the sum of s-channel resonances produces 
linear Regge trajectories in the t-channel [21]. This hypothesis is the basis 
of string models. 

The knowledge of Regge trajectories allows one to predict cross sections 
of two-body high energy processes with non-vacuum quantum numbers in 
the crossed channel. 

3.4 Regge pole theory for non- vacuum exchanges 

It has been suggested that reactions with non-vacuum quantum numbers in 
the crossed channel are dominated at high energies by the exchange of the 
Regge trajectories allowed by conservation of charge, spatial parity, G-parity 
and isotopic spin [22] ■ 

Assuming the dominance of the Regge pole contribution to the amplitude 
of diffractive processes, one obtains: 

da(a + b^a> + b>) = ^(t)^(t) /^ 2h( ' H) 
dt sm 2 (ira R (t)/2) \s 

+ ivJM^M (" x 2(afl,(t)_1) 

cos 2 (7ra R >(t)/2) \s 0/ 

+ 2 gaRa>(t)g bRb ,(t) g aR , a ,(t)g b R> b >(t) ( s^ <«*(t)+«*<tH») 

sm(ira R {t)/2) cos(vra iJ /( f) /2) \s 0/ 
x cos((7r/2)(a fl -ajy)). (36) 

This expression takes into account the exchanges of the trajectories with 
positive (R) and negative (R 1 ) signatures. Fitting Regge trajectories as a 
linear function of t gives a good description of the mass spectrum of the 
resonances belonging to the corresponding trajectories (see, e.g. Fig. [I]). The 
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assumption that the trajectories remain linear for t < leads to a reasonable 



description of the data using Eq. 36 



M 2 (GeV 2 ) 

7 
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a 4 (2040) >^/9 5 (2350) 



a 2 (1320) 
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p 3 (1690) 



Figure 1: The lowest mass mesons lying on the p Regge trajectory and on 
the nearly degenerate A 2 trajectory. 

The data exist for the cross sections of the following processes (see the 



review 



[16j 



Ti 



+ p — > 71° + n, 71 + p — >■ p(A 2 ) + n, K + p — > <f> + n, 



p + p — > A + N, etc. These data are described well at small t by the exchange 
of a few leading Regge meson trajectories. In the kinematics of large s and 
small u, the exchange of the baryon trajectories dominate such processes as, 
e.g., tx~ + p — >■ p + tx~ . The data are described well in terms of linear meson 
and baryon trajectories for —t < 0.5 GeV 2 (see, e.g., the review [16]). At 
the same time, for — t > 1 GeV 2 data can be interpreted as the evidence for 
flattening of meson and baryon trajectories to the values corresponding to 
the exchange of reggeized qq (qqq) systems [23j . Further experimental studies 
of high energy two-body reactions in the —t = const > 1 GeV 2 limit would 
be highly desirable. 



"Off w ^TFT^ ^T • (37) 



3.5 Non-universality of the Regge trajectories for the 
bound states containing heavy quarks 

The states containing heavy quarks belong to the Regge trajectories with the 

slope different from the one of the Regge trajectories for hadrons made of 

light quarks. We will give here two examples. 

(i) In the Mq — > oo limit, the masses of Qq states are Mgg^n) = Mg + A n , 

where A n are independent of Mq. This result follows from quarkonium 

models and is probably valid in QCD. Therefore, for linear trajectories one 

has: 

1 

2m q (a7- a ) 

This slope is different from the universal slope suggested in the string models 
for a hadron trajectory consisting of light quarks. 

(ii) For hadrons with hidden flavor, the slope of the QQ trajectory de- 
creases with the mass of the heavy quark in the Mq — > oo limit as: 

«io « ^ • ( 38 ) 

For an estimate, we used here quarkonium models with the Coulomb inter- 
action between quarks. 

To conclude, the large masses of heavy quarks supply extra scales in 
addition to the Aqcd scale, which suggests the existence of a variety of 
phenomena beyond the framework of the one-scale S- matrix theory [21] . For 
example, J /if) and T have significantly smaller radii than the pion: 

7V = 0.5 fm , rj/i/j fa 0.2 fm , r T ~ 0.1 fm . (39) 

As a consequence of color screening and asymptotic freedom, heavy quarkonia 
relatively weakly interact with hadrons made of light quarks. (The proba- 
bility of the pion field around a heavy quarkonium in the ground state is 
close to zero.) The observed total and partial widths and the cross section 
of diffractive photoproduction of J /if) at moderate energies are significantly 
smaller than those for vector mesons made of light quarks. 
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4 Pomeron theory of high energy soft QCD 
processes 

4.1 Introducing the concept of the Pomeron exchange 

We discuss briefly here the hypothesis of the Pomeron exchange dominance 
in the amplitudes of high energy processes. 
The Pomeranchuk theorem, 

<r(h + T) = (j(h + T), (40) 

was proven initially under the assumption that a to t{hN) —¥ const for s —¥ oo. 
The proof uses analytic properties of amplitudes in the energy plane and 
that the amplitude is predominantly imaginary at high energies. Indeed, 
A + (s,t) = c + [sln(—s) + u\n(—u)} ~ ic + ns which should be compared with 
A~ = c~[s ln(— s) — it In (—it) ] ~ 2c~sln(— s). The condition: A + ^> A~ 
requires that c~ = — the Pomeranchuk theorem [25] . In the more realistic 
case of a growing cross section, a weaker form of the theorem for s — > oo can 
be proven: 

a tot (/tH-T)Ax tot (/t + T)^l. (41) 

This theorem is confirmed by the data on c tot (pp) and <y tot (pp). 



As we explained in section 3J2 the behavior of the scattering amplitude 
A(s,t) = isf(t), which is typical for quantum mechanical problems with the 
absorptive interaction, contradicts the S-matrix unitarity relation for partial 
waves in the crossed t-channel. To resolve this contradiction, V. Gribov [15] 
suggested the behavior of the scattering amplitude (for large s and small t) 
of a general form that does not contradict the unitarity of the S'-matrix in 
the crossed channel: 

A(s,t) = is a{t) F(\n(s),t), (42) 

where F is a slow function of ln(s) rapidly decreasing with an increase of 
—t. The positive value of da(t)/dt at t = 0, which follows from the positive 
value of the partial waves (Im/;(s) > 0) and properties of the Legendre 
polynomials [15], leads to a decrease in the average —t for elastic scattering 
with an increase in s. 

Thus, the self-consistency of the theory requires that the radius of the 
hadron-hadron interaction should increase with energy. 

V. Gribov [2B] suggested that amplitudes of high energy processes are 
dominated by the special Regge pole trajectory exchange with the vacuum 
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quantum numbers in the crossed channel and calculated elastic and total 
cross sections. In parallel, G. Chew and S. Frautschi [22] also drew attention 
to the hypothetical Regge pole with a(0) = 1 and the vacuum quantum 
numbers that would be responsible for the forward scattering processes. This 
hypothesis reproduces the Pomeranchuk theorem and leads to an increase of 
the radius of the interaction with energy [TU [T5], I2S1 I2"7j : 

( s \ Qjp W / u \° ip ( i ) 
A(hT) = ghPh{t)9Tpr{t)- 



sin(7rajp(£)) 

9hiPh{t)9TPT [j-\ exp(i7rajp(t)/2) 
sin(7raj3(t)/2) 



(43) 



where ghPhif) and gxprif) are the residues of the Pomeron pole. The signa- 
ture factor follows from the symmetry of the amplitude due to the Pomeron 
exchange under the s -H- u transformation as explained in the previous sec- 
tion. 

4.1.1 Experimental evidence for the Pomeron trajectory 

The Pomeron trajectory is usually parametrized as a trajectory linear in t: 

a P (t) = ap + a'jpt . (44) 

The observed dependence of the total cross section of pp collisions on energy, 
which ranges from the fixed target energies to the highest energies currently 
measured at the LHC (■\fs=7 TeV), is well described by the Pomeron inter- 
cept [28]: 

a° P w 1.08 - 1.1 . (45) 

The generally accepted value of the slope of the Pomeron trajectory is 

a' P w 0.25 GeV" 2 . (46) 

The hypothesis of the dominance of the Pomeron exchange in high energy 
processes has found a number of experimental confirmations. Let us briefly 
outline basic discoveries. 

(i) The same Pomeron intercept describes the energy dependence of the 
total and elastic cross sections of pp collisions and of the total cross section 
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of exclusive photoproduction of p mesons off the proton target measured at 
HERA [28J . The cross section of exclusive photoproduction of J/ip mesons 
increases with energy more rapidly, see a recent summary in [29] (for the 
explanation, see the following subsections). 

(ii) The shrinking of the diffractive peak with an increase of the collision 
energy has been predicted in [221 I2E] and observed in pp and pp collisions (see 
the summary of the data in [3]]). The data for y/s < 1.8 TeV are consistent 
with the Pomeron trajectory being linear for —t < 0.5 GeV 2 . 

(hi) The first LHC data [31] report the t-slope corresponding to a faster 
rate of the shrinkage: for y/s between y/s = 2 TeV and \/s = 7 TeV, 
a'p(—t < 0.15 GeV 2 ) ~ 0.55 GeV -2 ; this value significantly increases with 
energy for larger —t. The interference between the single and multi-Pomeron 
exchanges produces a qualitatively similar behavior which we illustrate be- 
low. At small t, the amplitude for the single Pomeron exchange is conve- 
niently parametrized as: ImAi = cis ap ^~ l exp(Bt/2). Taking into account 
interference with the double Pomeron exchange, we obtain for the square of 
the ratio of the full amplitude to the single Pomeron exchange amplitude: 
|ImA/ImAi| 2 « 1 - 2(c 2 /ci)s 2a2p(4/4) ~ a2p(i) " 1 exp(- J Bt/4), where a are posi- 
tive and C2/C1 can be evaluated in the color fluctuation model. This formula 
produces the minimum which moves to smaller t with an increase of energies. 
This phenomenon has been discussed in the eikonal models, see, e.g., [32J; 
the problems of this approximation will be discussed in section [6] 

(iv) The global analysis of the world data on exclusive p-meson photopro- 
duction [30J gives a'jp = 0.126 ± 0.013(stat.) ± 0.012(syst.) GeV" 2 indicating 
non-universality of the t dependence of the effective Pomeron trajectory, see 
Fig. [2] The non-universality can also be seen from the comparison of the t 
dependence of ap(t) with the linear Pomeron trajectory that describes the 
pp data for the same energy interval. However the data do not exclude a 
possibility that a'jp for the p case is the same as for the pp case for —t < 0.15 
GeV 2 . Note here that selection of small t enhances the contribution of periph- 
eral collisions and hence suppresses the effects of multi-Pomeron exchanges 
and of blackening of the interaction at small impact parameters discussed 
in the next subsections. The data also indicate that ap > 1 for the entire 
studied t-range (—t < 1 GeV 2 ). In the case of elastic photoproduction of 
J/ip, ap ~ 1.2 for forward scattering and ap ~ 1 for —t ~ rn 2 j/^- Such a 
behavior is natural for the pQCD regime where double logarithmic terms are 
relevant for the interrelation between the t and s dependences [33] . 

(v) It is well known from the analysis of the ladder diagrams relevant 
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Figure 2: The data on the t dependence of the Pomeron trajectory as ex- 
tracted from p-meson photoproduction with a linear fit to the data and the 
Donnachie-Landshoff parametrization of the pp data. 



for the total DIS cross section [31] that a'p decreases with the virtuality 
of the external probe due to the suppression of the Gribov diffusion in the 
impact parameter space. The observed value of ot w for J/ip photoproduction 



is definitely much smaller than that for pp scattering, see Eq. 104. However, 



the errors in the J/ip case (see Eq. 104 ) do not allow one to establish whether 



a'jp extracted from the global analysis of the data for the J/tp case is smaller 
than for the p case. 

(vi) One of important confirmations of the dominance of the Pomeron 
exchange is the observation of the triple Pomeron diffraction [351 ES] at 
FNAL [37] and at the LHC [35], see the discussion in the next subsection. 

The dominance of the Pomeron exchange in the total cross section al- 
lows for a relationship between the cross sections of diffractive and inelastic 
processes [39]. Thus, the total cross section of diffraction is unambiguously 
calculable in terms of the cross section of inelastic processes, i.e., diffraction 
in high energy processes is a shadow of inelastic processes. 
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4.2 Pomeron calculus 

In this section we will discuss the physics of interacting Pomerons. The ideas 
and methods discussed in this section are now widely used in the evaluation 
of amplitudes of high energy processes, for a detailed discussion of the subject 
and proper references, see |17j . 

A single Pomeron trajectory generates multi-Pomeron branch points and 
related cuts in the angular momentum plane. This result follows from unitar- 
ity of the S-matrix for partial amplitudes continued into the angular momen- 
tum plane in the crossed t-channel [121 HO]- An exchange of n Pomerons leads 
to a branch point in the angular momentum plane located at [4"Tl IT2~| 140] : 

j{t)ajp(t/n 2 ) -n+l« n[a F (0) - 1] + (a/ F /n)t . (47) 

Realization of the Pomeron in the form of ladder diagrams helped to 
develop the diagram technique [12], which allows one to calculate the con- 
tribution of any number of Pomeron cuts to the cross section and, hence, 
helped to develop the Pomeron calculus. 

The simplest diagrams of the Pomeron calculus, which include a single 
Pomeron exchange, a two-Pomeron exchange, and a triple Pomeron exchange, 
are shown in Fig. [3j The coupling of two Pomerons to a hadron is expressed 
through the diffractive a + b — > X + b cross section that includes both the 
elastic and inelastic contributions. For t away from zero, the latter is en- 
hanced as compared to the elastic contribution since the t dependence of 
diffraction is weaker than that of the elastic scattering. This is because in 
quantum mechanics form factor of bound state decreases significantly more 
rapidly with an increase in -t than does the form factor describing sum of 
inelastic transitions. 




Figure 3: The single Pomeron exchange (a), the two-Pomeron exchange (b) 
and the triple Pomeron (c) diagrams. 
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If the intercept of the Pomeron were equal to unity (ap(0) = 1), all branch 
points would be located at j — 1 for t — 0. In this case, the calculation of the 
energy dependence of amplitudes of high energy processes leads to the scaling 
behavior of the Green's functions in the angular momentum plane in the 
vacuum channel. This technique developed in [43J helped to build the theory 
of second-order phase transitions and to analyze long-range fluctuations near 
the critical point. 

There exists a variety of experimental evidences for the important role of 
multi-Pomeron interactions: 
(i) The triple Pomeron diffraction gives the pattern of how multi-Pomeron 



interactions arise, see the discussion in subsection 4.4 



(ii) The shift to smaller —t of the position of the minimum of the elastic cross 
section with increase of energy indicates that the increase of the role of the 
multi-Pomeron exchanges with increase of s and t. 

(hi) The phenomenon of nuclear shadowing in hadron-nucleus collisions, 
where the incident hadron interacts with several nucleons of the nucleus, 
arises mostly due to multi-Pomeron exchanges. 



a *— I 1-» 

b > I L 



Figure 4: Ladder diagrams building the Pomeron exchange. 



4.3 Gribov diffusion in the impact parameter space 
within the Pomeron ladder 

The dynamical interpretation of the Pomeron exchange in the multiperiph- 
eral model or, equivalently, in the parton model, allows us to visualize the 
Pomeron in the phase space and to evaluate diffusion to large impact parame- 
ters. The Pomeron is modeled as the parton ladder (Fig. |4l) where transverse 
momenta of produced partons are independent of energy, [kf) = k^ [4*41 145] . 
Experimental analyses of the transverse momenta distributions of produced 
pions and kaons show that \/k\ « 0.3 — 0.4 GeV/c. However, the majority of 
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observed mesons originate from decays of heavier hadrons (direct hadrons) 
which have significantly larger average transverse momenta y k$ ~ 0.5 — 0.6 
GeV/c already at moderate \/s = 10 — 20 GeV. The distances in rapidity 
Ay between partons adjacent in the ladder are independent of the collision 
energy. (The data on inelastic hadron production in high energy scattering 
find significant rapidity correlations between hadrons only for Ay < 1.) In a 
naive picture of the Pomeron described by a single ladder, each decay corre- 
sponds to a step of a random walk in the b space of the length oc 1/ko and 
the number of steps oc ln(s / /i 2 ) / Ay . This leads to an increase of b 2 with 
energy, b 2 oc ln(s//x 2 ), and hence to finite a' |Hj. 

Diffusion to large impact parameters manifests itself directly in elastic 
hadron-hadron collisions. In the impact parameter space, the imaginary 
part of the partial amplitude due to the Pomeron exchange has the following 
form: 

d 2 q t *A(s,<g) cs ap ^- l exp{-b 2 /2B) 



lmf{s, b 2 ) = / — -^ exp(iq t ■ b) 



2s B 



where B = B + 2a'jp ln(s/so) is the t-slope of the elastic cross section. We 
parametrize the product of the residues as cexp(B t/2) to take into account 
that Pomeron dynamics is dominated by peripheral, that is, small t processes. 



It follows from Eq. 48 that the impact parameters essential for elastic pp 



collisions increase with energy as: 

<6 2 > el ^ B + 2t4ln(s/so). (49) 

In the case of the total cross section, 

<" 2 >« = 2<"V (50) 

Thus the assumption of the dominance of a single Pomeron exchange 
implies the dominance of peripheral collisions in hadron-hadron interactions 
at ultra high energies. 

If spin-flip effects are neglected, the amplitude A(s,t) can be extracted 
from the differential cross section within the forward peak: 

da el /dt=— \A(s,t)/s\ 2 . (51) 

l07T 

In elastic pp collisions, spin-flip effects are small for — t < 1 GeV 2 [IB]. There 
exists an insignificant correction due to the possible ambiguity in the sign 
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of the amplitude beyond the cross section minimum, which for yfs = 7 TeV 
occurs at —t ~ 0.5 GeV 2 . Thus the impact parameter distribution can be 
measured using the elastic cross section data. 

The parton ladder model of the Pomeron allows us to evaluate properties 
of multi-Pomeron exchanges and Pomeron loops. To explain the role of multi- 
Pomeron exchanges, let us consider a simple model where the product of the 
Pomeron residues is parametrized as above by the factor of gh 1 Ph 1 gh 2 iPh2 ex P B/2t. 
In this case the impact parameters corresponding to an exchange of n Pomerons 
rapidly decrease with n: 

(#>«<&>/n, (52) 

where (bjp) = B are the impact parameters characteristic for the single 
Pomeron exchange; B is the slope of t dependence of the elastic cross section. 
(The inclusion of inelastic intermediate states for multi-Pomeron exchanges 
would result in a slower t dependence and, hence, in the further reduction of 
the average b 2 for multi-Pomeron exchanges.) 

The relative contribution of the multi-Pomeron exchanges as well as the 
effects of the Pomeron self-interactions grow with an increase of energy. As 
a result, the amplitude of the elastic pp collision becomes completely absorp- 
tive, Imf(b, s) «1 for small impact parameters, see the discussion in the 
next subsection. The blackness of the interaction at small impact parame- 
ters somewhat suppresses diffusion in the impact parameter space because 
the trajectories entering the absorption region disappear from Gribov diffu- 
sion. 

Another effect leading to the suppression of Gribov diffusion is an increase 
with collision energy of the probability of hard processes with the correspond- 
ing tendency of the disappearance of significant differences between soft and 
hard QCD processes. It is rather difficult to observe this phenomenon di- 
rectly since the majority of observed pions results from the decays of heavier 



resonances (section 4.3) which have average transverse momenta as high as 



0.6 GeV/c already at fixed target energies. 

4.4 Observation of multi-Pomeron interactions: soft 
diffraction in the triple Pomeron limit 

Probability of the processes with a large gap in rapidity can be evaluated 
using the Pomeron exchange. The Pomeron calculus predicts diffraction of 
an incident hadron into hadronic states whose invariant masses Mx are large 
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and the rapidity gap is large and increasing with energy, 

h + T -> M x + rapidity gap + T' , (53) 

in the kinematics where M x /s = Xp is small and constant and s — > oo. 
Here 1 — x p is the fraction of the target momentum carried by the final state 
hadron T". Since Mx is large and increases with s, the sum over diffractively 
produced states can be substituted by the Pomeron exchange for xp <C 1. 
Thus the process of diffraction of an incident hadron into a large mass state 
probes the triple Pomeron vertex (the process corresponding to the diagram 
in Fig. 3c). 

The observation of this process at FNAL [37] and the LHC [38] is a direct 
demonstration of how the Pomeron branch points arise and of their important 
role in high energy processes. 

The consideration of the diagram c in Fig. [3] allows one to predict the 
dependence of the cross section on xp 



da{h + T -> "M x " + rapidity gap + T')/dx P d 2 q t = (54) 



2 



where q t is the transverse momentum transferred to the target. A distinctive 
feature of this formula is the singularity of the cross section at qt = and 
xp —¥ 0, if the triple Pomeron vertex g%p{qt = 0) is different from zero. If 
cvip(0) > 1, this singularity will be present at non-zero t as well. This sin- 
gularity lies in the unphysical region since xp = requires infinite energies. 
With an increase of s, multi-Pomeron exchanges and Pomeron loops, which 
were effectively forbidden at lower energies by energy-momentum conserva- 
tion, become progressively important. 

Inelastic diffraction can occur only at large impact parameters since in- 
elastic processes fill the rapidity gap at small impact parameters, where the 
interaction is practically completely absorptive. The Pomeron calculus takes 
this effect into account by including multi-Pomeron exchanges that strongly 
screen the triple Pomeron contribution, especially at small impact param- 
eters. The overall effects are the reduction in the large-mass diffraction, a 
faster decrease in the cross section with gf, and a reduction of a e p{t ~ 0). 
Screening of the triple Pomeron vertex was evaluated in the generalized 
eikonal approximation in [37] and was found to be very large. The behav- 
ior of inelastic diffraction in pp collisions that we described above has been 
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observed at FNAL [37]. Most of the FNAL data [37] correspond to rela- 
tively large xp > 0.01 so that the contributions of secondary trajectories 
play an important role. The recent LHC data [38j observed similar regulari- 
ties but at smaller xjp, where the contribution of secondary trajectories can 
be neglected, while the screening effects may play a role; a^ ~ 1.05 was 
reported [35] . 

4.5 Blackening of hadron— hadron interactions at cen- 
tral impact parameters 

Probability conservation, i.e., unitarity of the S matrix in the s-channel, 
restricts the high energy behavior of the total cross sections of hadronic 
collisions and cross sections of diffractive processes: 

1 2 

ha f(b,s) = - |/(fo, s) \ + positive terms. (55) 

High energy processes are predominantly inelastic so that the partial waves 



are predominantly imaginary. Indeed, it follows from Eq. 55 together with 



the expression for the cross section of inelastic processes [48J that: 

&mei(s,b) = l-\S l - 1| 2 , (56) 

where Si = 1 + if] is the matrix element of the S'-matrix corresponding to 
the orbital momentum I. If Im/(6, s) exceeds unity, a- me \(s, b) would start 
decreasing with an increase of Im f(b, s) in contradiction with the dominance 
of inelastic processes. As a result, one concludes that at high energies the 
partial waves for elastic collisions cannot exceed unity: 

Im/(6,s)<l. (57) 

The upper boundary for the total cross section follows from the S-matrix 
theory [7] . Indeed, if in the kinematical region of t restricted by the singu- 
larities of the amplitude in t plane the scattering amplitude increases with 
energy not faster than a polynomial, i.e., if ImA(s,t) < s N , then 

a tot <c\n 2 (s/s ). (58) 

The formal derivation [49] uses analytic properties of the amplitude in the 
t-plane and unitarity of the S'-matrix to derive the polynomial boundary on 
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the amplitude that we mentioned above. This derivation also allows one to 



evaluate the coefficient c in Eq. [58] whose value turns out to be unrealistically 
large, being significantly larger than that arising from fits to the pp data. 
As we explain above, the amplitude A(s, t) can be unambiguously ex- 



tracted from the data on elastic pp collisions using Eq. 51 In the kinematics 
achieved at FNAL for pp scattering, cr to t{pp) ~ 167rI?.This means that for the 
pp interaction at the zero impact parameters, the partial amplitude is close 
to unity. If we define (for the illustration purposes) proximity to the black 
regime as a condition that probability of the inelastic interaction is > 0.75 



( see Eq. 56), the interaction will be close to being black for b < 1.0 fm at 
y/s = 7 TeV and for ~ 15% smaller values of b at yfs = 2 TeV. For larger 
b, which dominate in the inelastic cross section (the median b ~ 1.4 fm at 
-y/s = 7 TeV), the interaction is grey and rather far from the black regime. 

It is well known from textbooks that the assumption that the interaction 
is black at all impact parameters leads to the following two predictions: 

<j{hT) el /a{hT) tot = l - (59) 

and 

da/dt{pp — > diffractive state + rap. gap + p)t=o = . (60) 

The second property of the black disc regime (BDR) follows from orthogonal- 
ity of the wave functions of the eigenstates of the Hamiltonian corresponding 
to different eigenvalues. 



Eq. 59 is in the evident disagreement with the data as the u c i/a tot ratio at 
the LHC, while slowly increasing with energy, still reaches only the value of 
about 25%. Also, at the LHC the ratio of inelastic and elastic diffraction is 



close to unity in the apparent contradiction with Eq. 60 Taken together, the 
discussed features of the data point out that, effectively, the pp scattering at 
FNAL and the LHC has a small black spot at small b surrounded by a large 
grey area. The grey part mostly represents the contribution of peripheral 
collisions to the total cross section and it is dominated by exchanges of the 
Pomeron and Pomeron cuts. 
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5 Space— time evolution of high energy pro- 
cesses 

5.1 Introduction 

The phenomenon of the linear increase of longitudinal distances with an in- 
crease of the collision energy was first understood in QED for the propagation 
of charged particles through the medium and led to the so-called Landau- 
Pomeranchuk-Migdal effect [50]. In QCD, the linear increase of longitudinal 
distances with the collision energy [51] is a fundamental property of high en- 
ergy processes. This and following sections discuss new coherent phenomena 
that follow from taking into account the space-time evolution of high energy 
processes in QCD. Note that the concept of the space-time evolution of high 
energy processes [52] is beyond the framework of the S'-matrix approach. The 
aim of this section is to explain some properties of the space-time evolution 
of the scattering processes and their implications. 

5.2 Linear increase of longitudinal distances in high 
energy processes with energy 

It follows from the application of the energy-time uncertainty principle to 
scattering processes in QCD that when the energy of the projectile in the tar- 
get rest frame, E, is large enough, the quark-gluon configurations satisfying 
the condition: 

z « (1/(A£)) « 2E/(M 2 - m 2 h ) > R T (61) 

are formed before the target and these configurations are frozen before the 



collision. In Eq. [6TJ (M 2 - m 2 hadron ) = ^(m 2 . + &&)/* ~ m Lw) and z i 
is the fraction of the incident particle momentum carried by the constituent 
i. Thus, the space-time evolution of the scattering processes differs from 



that in non-relativistic physics. Eq. 61 directly follows from the Lorentz 



transformation applied to projectiles with an energy-independent number 
of constituents. One can easily check the validity of the above estimate 
of longitudinal distances by analyzing the dominant ladder diagrams using 
popular approaches to high energy processes. 

Here we will present reasonings applicable also in the non-perturbative 
QCD regime. Let us consider the virtual photon-target scattering in the 
deep inelastic (Bjorken) limit: — q 2 — > oo, —q 2 /2(q ■ px) = const, where q 
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and pt are the four-momenta of the photon and the target, respectively. It 
follows from the optical theorem that the total 7*-target cross section has 
the following form: 

a = -Iim4( 7 *+T^ 7 *+T) = - f exp{iq-y) (T\ [J,(y), J M (0)] \T) d*y . (62) 
s s j 

As a consequence of causality, only the y 2 = t 2 — z 2 — y 2 > region contributes 



to Eq. 62 At large energies, go = yH ~ Q 2 ) ~ Qz — Q 2 /2q z and, thus, 
i(q ■ y) « i(qo(t — z) — zQ 2 /2q ). Since the Q 2 dependence is contained only 
in the second term and the cross section decreases with an increase of Q 2 , 
the direct analysis of the representation of the cross section in the form of 



Eq. 62 shows [5T] that essential distances in the integral are t « z and 

z « 2q /Q 2 . (63) 

5.3 Cancellation of the contribution of planar/ Glauber- 
approximation diagrams 

It has been understood long ago that the Glauber (eikonal) approximation — 
being a very popular method of modeling of high energy processes in nuclear 
and particle physics — is actually inapplicable beyond the non-relativistic do- 
main. The dominance of large longitudinal distances changes qualitatively 
the pattern of multiple interactions. 

In the non-relativistic quantum mechanics, the eikonal approximation fol- 
lows from the Schrodinger equation when the kinetic energy of the incident 
particle significantly exceeds the potential of the interaction J3S]. In the 
Glauber approximation, high energy interactions of the projectile with a tar- 
get occur via consecutive rescatterings of the projectile off the constituents 
of the target. The projectile is on its mass shell between the interactions — 
one takes the residue in the propagator of the projectile (Fig. [5k). However, 
the Glauber approximation contradicts the QCD-based space-time evolution 
of high energy processes dominated by particle production. Indeed, as the 
essential distances become significantly larger than the distances (time in- 
tervals) between consequent rescatterings [51], there is no time for a frozen 
configuration in the projectile to recombine into the projectile during the 
time of the order of Rt since it is much shorter than the lifetime of the 



configuration (Eq. 61). 
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Figure 5: (a) The planar diagram for double scattering, (b) A non-planar 
diagram for double scattering. 



There are two independent theoretical proofs that the contribution of pla- 
nar diagrams to the double Pomeron scattering amplitude is actually zero. 
It was found that in the case of high energy scattering in a quantum field 
theory, the contribution of the planar diagrams with intermediate states cor- 
responding to the projectile on its mass shell drops with the incident energy 
as 1/s. Indeed, the integral over the square of the mass of hadrons pro- 
duced in the iP-hadron collision, M 2 , is zero in the case when the Feynman 
diagrams have only s or u cuts because the contour of integration can be 
moved in the direction where there are no singularities in M 2 [1TI W2\ . The 
integral over the large circle is zero as a consequence of a decrease of the 
amplitude with M 2 . The eikonal diagram (Fig. pk) belongs to the class 
of Feynman diagrams where the cancellation occurs. Crossed (non-planar) 
diagrams (Fig. [5b), which have cuts both in s and u, give a non-zero contri- 
bution. 

Taking in account energy-momentum conservation leads to the same con- 
clusion [53] . Indeed, the eikonal diagrams correspond to an inelastic interme- 
diate state described by the Pomeron exchange at the double energy 2s. On 
the contrary, in the crossed diagram the energy is divided between constituents 
before the collision. If one parton carries the fraction z of the incident hadron 
momentum and another parton carries the fraction z' < 1 — z, the total en- 
ergy of the produced hadronic state is sz + sz' < s. Both arguments can be 
easily generalized to the case when the wave function of the initial hadron 
contains many constituents. 

Using the technique of the Pomeron calculus, V. Gribov showed that in a 
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quantum field theory, the contribution of the non-planar diagrams relevant 
for the multiple Pomeron exchanges to the total cross section can be rewritten 
as a sum of the eikonal term and the inelastic diffraction contribution. The 
resulting Gribov-Glauber model [54J is in agreement with the data on nuclear 



shadowing in hadron-nucleus interactions [55]. We will explain in section 6.4 
how taking into account color fluctuations allows one to evaluate the relative 
contributions of multiple scatterings in the Gribov-Glauber model. 

6 Fluctuations of color in diffractive phenom- 
ena 

6.1 Introduction 

QCD predicts new types of diffractive phenomena as compared to the S- 
matrix approach since the wave function of an energetic incident hadron is 
formed long before the target and the transitions between different config- 
urations in the wave function occur at distances comparable with the char- 
acteristic longitudinal distance called the coherence length. Therefore, the 
cross section is calculable in terms of the instant quark-gluon configurations 
in the projectile (see the discussion in the previous section). In the exclusive 
processes, where the incident hadron is squeezed in the transverse direction 
by the choice of the specific final state, the spatially small wave packet of 
quarks and gluons weakly interacts with a target in a rather large interval 
of collision energies. This phenomenon is calculable for hard diffractive pro- 
cesses in the form of the special QCD factorization theorem [Ml 155"] and has 



been observed in a variety of experiments (see section 7.7 and the review 
and references in [57]). Thus QCD predicts fluctuations in the strength of 
the interaction since the interaction differs for different configurations of con- 
stituents in the wave function of the incident hadron. This is an intuitive 
justification of the necessity to use the concept of the distribution over cross 
sections instead of the average cross section. 

6.2 Suppression of the strong interaction due to screen- 
ing of color 

At large energies the wave function of the incident hadron is formed before 
the target and is frozen, if (2E / AM 2 ) 3> Rt- Different configurations of 
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constituents in the wave function of the incident hadron interact with the 
target at different strengths. This is an important property of QCD where 
the interaction is proportional to the area occupied by color since the color 
charge of a hadron is zero. For illustration purposes we begin with quark 
models of hadrons and then derive formulas in QCD. 

A popular model was suggested by Low and Nussinov [55] 159] in which 
the total hadron-hadron cross section is described by the exchange of two 
gluons. Low further argued that the cross section is proportional to the 
region occupied by color in the hadrons: 

<r{hT) = crf, (64) 

where r t is the transverse radius of the smaller hadron. The derivation in- 
volves taking into account the gauge invariance, the zero color charge of 
hadrons as well as an implicit assumption that the momenta of constituents 
within the hadrons are significantly larger than the transverse momenta of 
the exchanged gluons. Eq. [64] was elaborated on in the constituent quark 



model with a two-gluon exchange between h and T [BO]. Eq. 64 can be ques- 
tioned since in the model, the average size of configurations involved in the 
scattering is comparable to the scale of non-perturbative QCD phenomena 
and the restriction by a two-gluon exchange cannot be justified. 



Eq. [64] can be reformulated to include the full QCD. If the incident me- 
son is in a quark-gluon configuration whose transverse size is significantly 
smaller than the scale of non-perturbative QCD phenomena, ^Aq CD <C 1, 
the application of the technology of the QCD factorization theorem [56J and 
the QCD evolution equation for parton densities allows one to calculate the 
cross section of the hadron interaction with a target T. The derived expres- 



sion [6 1 » 62 , 63J also contains the factor of r\ as in Eq. 64 which represents the 
coordinate space equivalent of approximate Bjorken scaling for DIS processes. 
However, in addition, the final expression contains the factor xGt — the gluon 



distribution in the target T absent in Eq. [64 



..FV 



= V h {n = 0,r z = 0)—r 2 t a s (Q 2 )xG T (x,Q 2 = — 2 ) . (65) 

Here F 2 is the Casimir operator for the quark (gluon) dipole; r t is the trans- 
verse distance between the quark and the antiquark; iph{ r t) is the small trans- 
verse size component of the incident hadron wave function; x ~ l/(4srf ). 
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Eq. 65 can be obtained also from the formulae derived in [51] in the leading 
a s ]xl(xq/x) approximation. In the case of nucleon one needs to take into 
account that the square of nucleon wave function at small ru is proportional 
to r| and to make the substitution 

r 2 -)• (r*i - (r t2 + r t3 )/2) 2 + (r t2 - (r a + r t3 )/2) 2 + (r t3 - (r tl + r t2 )/2) 2 . (66) 

In this case, the three quarks act as a symmetrized superposition of three 
dipoles stretched between one quark and the transverse center of mass of the 
other two quarks. 

6.3 Perturbative Pomeron 

The asymptotic behavior of the amplitudes of high energy processes in the 
vacuum channel in a quantum field theory containing vector particles was first 
investigated in QED. In [H3] the sum of the leading a 2 em ln(s/s ) terms in the 
cross section of the photon-photon scattering was calculated. (Note that in 
QED, in the lowest-order over coupling constant, this cross section is inde- 
pendent of energy.) This idea was applied to QCD in j66l|67] to study ampli- 
tudes of high energy processes in the kinematics where ln(x /x) ^> ln(<5 2 /Qo) 
by summing leading a s ln(xo/x) terms — the so-called perturbative Pomeron. 
A priori one can try to justify this approximation in the case of scattering of 
two small dipoles of the transverse size oc 1/Q within the restricted kinemat- 
ical domain of rapidities y < yo(Q 2 ,x), where the coupling to the colliding 
dipoles is perturbative. The kinematical boundary for the applicability of 
this approximation — yo(Q 2 , x) — as well as for the decomposition over powers 
of 1/Q 2 arises due to diffusion in the space of transverse momenta to the 
non-perturbative domain [68] . 

In the leading log approximation, the cross section grows as: 

cTdipoic-dipoic oc {l/xf , (67) 

where 

^ ^■ 4 " 12 . (68) 

For Q 2 ~ 2 GeV 2 , a s ~ 0.25 leading to /3 « 0.7. Note that the actual formula 



derived in [66J, EZ] is significantly more complicated than Eq. 68, which is just 
a popular fit to this formula. The derived expression [OH] corresponds to a 
cut in the angular momentum plain so that it involves a mathematical object 
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that is different from the Pomeron Regge trajectory discussed in section |4j 
It becomes a sum of the poles in the angular momentum plane in the large 
N c approximation. 

Assuming that the number of radiated gluons is sufficiently large, the 
diffusion equation was derived for the motion in the plane of ln(p 2 /p 2 ). It 
was observed that diffusion both to large and small p t is present [66, 69J. 
Significant diffusion to small pt, i.e., into the non-perturbative domain raises 
questions about the validity of pQCD approach because the answer depends 
on the treatment of the badly understood infrared region. The main dif- 
ference between the perturbative Pomeron and the Pomeron trajectory dis- 
cussed in section El is significant diffusion to large parton momenta (in addi- 
tion to diffusion to small parton momenta). Such diffusion is absent in the 
non-perturbative Pomeron which is modeled by the non-perturbative ladder 
discussed in section HI 



The NLO correction to Eq. 68 was found to be so large [701 [71] that it 
dominates the LO expression for a wide range of a s , which leads to (5 ~ 0. 
This is primarily because the LO as well as NLO approximations ignore 
energy-momentum conservation. The poor convergence of the series for the 
total cross section in terms of powers of ln(s//i 2 ) was first demonstrated in 
QED [72] by the direct calculation of the lowest order diagrams for the e + e~ 
pair production in electron-electron scattering. For the contribution of the 
dominant two-photon mechanism, it was found that 

a = a^ m c(1.041n 3 (s//i 2 ) - 6.591n 2 (s//i 2 ) - 11.81n(s//i 2 ) + 104 + 0(fx 2 /s)) , 

(69) 
where /i is the electron mass. It was explained in [72] that a fast growth of 
the coefficients in front of the powers of ln(s/so) reflects the highly restricted 
phase space for obtaining logarithmic contributions. 

The resummation models [731 El] more smoothly match with the formulas 
of the DGLAP approximation where the conservation of the longitudinal 
component of the momentum is exact. 

The same approach is often applied to scattering of a small dipole of the 
diameter d rj 1/Q off the nucleon whose diameter is ~ l/{2m 7T ) 1 which is 
relevant for such processes as inclusive DIS, exclusive production of vector 



mesons, etc. Although the energy behavior given by Eq. [67] obviously con- 
tradicts the data, resummation approaches can fit the data since they lead 
to the results close to those obtained in the DGLAP approximation. Up to 
now the resummation approaches have not been applied to the description 
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of hard diffractive processes. 



One should note that the probability conservation in the form of Eq. |118 



is also violated in the resummation models at sufficiently small x and small 
impact parameters (see the discussion in section |8J). 



6.4 Distributions over strengths of the interaction for 
hadron and photon projectiles 

Before QCD was recognized as the theory of the strong interactions, in the 
framework of the parton model in which the strength of interactions is pro- 
portional to a number of wee partons in the projectile configuration, Pumplin 
and Miettinen suggested the description of high energy diffractive processes 
in terms of the probability distribution over cross sections, Ph{cr) [ZS]. It was 
understood later on that such a distribution originates from the dependence 
of the cross section on the instant transverse radius of the color distribution 
in the incident hadron. 

Constructively, -P/i(cr) is defined in terms of its moments, 

(a k )= [ P h (a)a k da, (70) 



with additional general QCD restrictions on the form of Ph{cr). The case of 
k = corresponds to the normalization condition for Ph{<j). By definition, 
the first moment of Ph{&) (k = 1) is the total cross section of hN scattering. 
The dispersion of the distribution over a is given by the ratio of inelastic and 
elastic diffraction at t = [75J: 

d<r(h + N^h' + N)/dt\ t= o (a 2 ) - (a) 2 _ 
da(h + N^h + N) ]t=0 "" (a) 2 " ^ ' l ' 

where the state h! differs from the state h. 

The behavior of Phi?) at small o follows from the kind of quark counting 
rules. Taking into account the number of valence quarks in a hadron h and 
using approximate Bjorken scaling in the form explained above (a ex rf), it 
was found that [6H [76] : 

P h {a -)■ 0) oc a K+"?+^- 2 ) , (72) 

where rii is the number of valence constituents in the incident hadron in the 
configuration participating in the scattering proces^j The data on diffractive 

1 In the analysis of |75] it was assumed that Pn{o) contains the term oc 6(a). 
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p+ 2 H — > X+ 2 H scattering provides an additional constraint, ((cr — (cr)) 3 ) ~ 
at t/snn = 30 GeV [76J. The information on the first three moments of 
P/j(cr) and its behavior at a — >• allows one to reconstruct the form of Ph(&) 
for the pion and nucleon projectiles (Fig. [6]). 




Figure 6: The distribution over a for fixed target energies extracted in [76J. 
The set of curves reflects uncertainties in the extraction procedure. The 
rectangular area in the upper left corner is the pQCD evaluation of P n (cr — > 

At present, the distributions are reconstructed for the energies of incident 
hadrons in the range of a few hundred GeV in the target rest frame. With 
an increase of the energy, the edge of the distribution moves to the right 
diminishing probability of weakly interacting configurations. The dispersion 
of the distribution grows with y/s up to \/s = 50 — 100 GeV, where u a ~ 0.3 is 
reached, and it starts dropping for larger \fs. Preliminary LHC data indicate 
that for pp scattering at \fs=7 TeV, u ff ~ 0.2. Note that uj a ~ 0.2 — 0.3 
corresponds to very large fluctuations of the strength of the interaction. For 
example, if one models P/v(a) as a superposition of two scattering states, cr's 
for these states would be a to t(l ± \/^7) corresponding to o\ ~ 55 mb and 
o"2 ~ 145 mb at \fs = 7 TeV. 

In the case of photon, 



P 7 (a-^0) oc l/a 



(73) 



which follows from the presence of point-like qq configurations in the photon 
the wave function. 
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The Gribov-Glauber model allows for a transparent interpretation in the 
formalism of cross section eigenstates of Good and Walker [77] . Indeed, 
different configurations are absorbed with the strength given by the Glauber 
model for a given a and incoherently contribute to the total cross section [78J . 

The concept of Ph{&) allows us to build a compact implementation of the 
Gribov-Glauber series for the total cross section of hA scattering: 

a tot (hA) = J daP h (a) J d 2 b2[l - e "^( 6 )/ 2 )] . (74) 

It allows us also to calculate the total cross section of inelastic coherent 
diffraction off nuclei (hA — > h'A) in a good agreement with the data J79J EQl 
I8T] . It also allows us to model deviations from the Glauber model in inelastic 
proton (nucleus)-nucleus collisions [82]. Note here that deviations from the 
eikonal approximation for the interactions with j nucleons is given by the 
(<r J ) / (a) J ratio, which rapidly grows with j. 

Knowledge of the first three moments of the distribution over a is suf- 
ficient to describe many nuclear phenomena, in particular, the total and 
coherent inelastic diffraction cross sections. The evaluation of more compli- 
cated phenomena such as, the tail of the hadron multiplicity, at present is 
model-dependent since the behavior of -P/i(cr) at large a is far from being 
understood. 



6.5 Diffraction in deep inelastic collisions as a pattern 
for the fluctuations of color 

In this subsection we will consider diffraction in deep inelastic ep scattering: 

7* + p — > X +rapidity gap +p (75) 

in the Bjorken limit. 

The standard picture of DIS is that of the absorption of the virtual photon 
by a parton (quark or antiquark) that carries a fraction x of the light-cone 
momentum of the nucleon with radiation of gluons in the initial and final 
states. In such a picture, which was quite popular before the first measure- 
ments at HERA, the hard gluon radiation should fill the whole available 
rapidity interval (in addition, partons should be emitted to screen the delo- 
calized color in the final state) and lead to the disappearance of diffractive 
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processes. However, such processes were observed with a significant proba- 
bility even at very large Q 2 . 

We will focus our attention on the limit when M\/Q 2 = const. In this 
limit, it is convenient to introduce the variable j3: 



Q 2 + M 2 X 



?=n^- ( 76 ) 



The variable (3 is related to the fraction of the momentum lost by the nucleon, 
xjp, as 

/3 = x/ X]P . (77) 

It is convenient to introduce " conditional" or fracture structure functions 
for the processes where one hadron is fixed in the fragmentation region. For 
diffractive processes, one usually uses the notation F p '(/3, Q 2 , xp, t) and one 
can also introduce quark and gluon diffractive parton distribution functions 
(PDFs) that depend on the same variables. Since hard processes occur locally 
in transverse momenta and rapidity, the increase in the resolution should not 
affect the properties of the nucleon fragmentation region. Hence one should 
expect that the diffractive PDFs should satisfy the same DGLAP evolution 
equations as usual PDFs. 

Extensive data on various large-mass diffractive processes have been ob- 
tained at HERA (for the recent results, see [HS]). The principal findings are 
the following: 

• The leading twist approximation with the same diffractive PDFs consis- 
tently describes the Q 2 evolution of the inclusive diffractive cross section and 
the diffractive cross sections of dijet (X= jeti + jet2 +X') and heavy flavor 
production for fixed xp. Factorization was formally proven in [84"] . 

• The data are also consistent with the Pomeron factorization: 

ff D (j3,Q 2 ,x P ,t)=r(xp,t)f j (x,Q 2 ). (78) 

The xp dependence of r(xp,t) is given by the same expression as that for 
soft diffraction, which employes ap(0) = 1.11, which is close to op(0) ex- 
tracted from the analysis of soft diffractive processes, total pp cross sections, 



and exclusive light vector meson photoproduction (subsection 5.2). The ob- 
served value of the intercept is significantly smaller than ap(0) for such hard 
exclusive diffractive processes as J /if) photo(electro)production. It was also 
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(79) 



found that gluons play a very important role in diffractive dynamic: 

nu 4 

£«* Jo #/*/«(& Q a ) 

for Q 2 ~ a few GeV 2 . 

• The overall probability of diffraction in DIS 

#(x, Q 2 ) = a diff (x, Q 2 )/a D i S (a;, Q 2 ) , (80) 

is of the order of 10% and grows with a decrease of x for fixed Q 2 . 

As we mentioned above, hard processes cannot screen the quark (anti- 
quark) emitted by a highly virtual photon. Therefore, pQCD states should 
contain no rapidity gaps and, thus, diffraction should be part of the non- 
perturbative initial condition for the QCD evolution equation, which is far 
from trivial to implement in the infinite momentum frame. The presence of the 
leading twist diffraction imposes constraints on sea quark and gluon nucleon 
PDFs at the starting point of the evolution — they should exhibit the small x 
behavior consistent with the soft Pomeron limit, i.e., they are allowed to grow 
only slowly with a decrease of x. 

The significant value of the cross section of diffraction can be understood 
in the formulation of the parton model in the target rest frame suggested 
by Bjorken in 1970 twenty years before diffraction in DIS was observed. At 
small x in this reference frame, the virtual photon transforms into qq pairs 
well before the target (section [6]) . To satisfy Bjorken scaling, it is necessary 
to assume that only the qq pairs with k t < k t o and the light-cone fractions 
satisfying the condition k 2 /z(l — z) ~ Q 2 should interact with the target 
at small x with a strength comparable to that of the pion-nucleon interac- 
tion, while the contribution of the pairs with k t 3> k t0 should be strongly 
suppressed. The low k t pairs are aligned along the photon direction — hence 
they are referred to as the aligned jet model (AJM) [85]. The 1/Q 2 behavior 
of this contribution is due to the small phase space allowed for these config- 
urations. In the coordinate space this corresponds to production of a qq pair 
at a distance 2q /Q 2 from the target and the expansion of the pair to the 
hadronic-scale size of l/k t0 by the time it reaches the target [86j. 

In QCD, the parton picture is modified by the following two effects. 
First, the emission of a large-size qq pair without the associated gluon emis- 
sion is suppressed by the Sudakov form factor. Inclusion of this emission 
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leads to the scaling violation, but it does not change the size of the quark- 
gluon configuration — the QCD AJM [87]. Second, while the interaction of 
qq pairs with the large transverse momenta up to k t oc Q is suppressed 
by the 1/Q 2 factor due to color transparency, it also contains the factor 
a s (Q 2 )xGN(x,Q 2 ). As a result, there is conspiracy between the hard and 
soft contributions — both of them are oc 1/Q 2 , with the hard contribution 
being numerically suppressed at moderate x ~ 1CT 2 , but gradually grow- 
ing in importance with a decrease in x due to the corresponding increase 
in xGn- The contribution of large masses 3> Q 2 , i.e., |3 C 1 (the triple 
Pomeron processes) requires sufficiently small x to reveal itself. Otherwise, 
it is suppressed by energy-momentum conservation. 

Thus, the probability of diffraction in the aligned jet model (DGLAP 
approximation) is comparable to that in hadron-nucleon scattering. At the 
same time, the contribution of small-size configurations to the cross section 
of diffraction (integrated over (3) is suppressed relatively to the inclusive cross 
section by the factor 

^diff a 2 s (Q 2 )(xG T (x,Q 2 )) 2 /BQ 4 2 n^/nrfl (ai\ 

=c , /n 2 wrm oca s {Q )xG T {x,Q )BQ , (81) 

Ctot |hard a s xG T {x, Q 2 ) I Q 1 ) 

where B is the slope in the t dependence of the diffractive cross section. This 
ratio rapidly decreases with an increase in Q 2 and increases (for fixed Q 2 ) 
with a decrease in x. In the fast frame, AJM configurations are equivalent 
to the presence of local (in rapidity) color screening of qq pairs in the small 
x nucleon wave function [88] . 

The observed soft xp dependence of diffraction is natural in the QCD 
AJM in the kinematic range of the validity of the DGLAP approximation. 



The approximate soft factorization (Eq. 78 ) is natural in the ladder mod- 
els of the Pomeron since in these models, the structure of the ladder does 
not depend on the rapidity for op(0) ~ 1. An increase of the probability of 
hard small x processes would result in breaking of the soft factorization at 
very small Xp. 

At sufficiently small x and moderate Q 2 , the hard contribution may be- 
come significant. Attempts to incorporate these higher twist contributions 
were taken in a number of the dipole models of the j*N interactions (see, 
e.g., [89] and references therein). However, most of these models ignore the 
Q 2 evolution of the AJM component. 

To conclude, the significant cross section of diffraction in DIS is another 
demonstration of the important role of color fluctuations in the virtual pho- 
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ton wave function and of the dominance of soft Pomeron physics even in 
seemingly hard processes. In such processes soft dynamics together with Q 2 
evolution gives a significant contribution to the total cross section. 

7 Hard Exclusive Processes 

7.1 Introduction 

It has been understood since nearly two decades ago that a number of two- 
body and quasi two-body processes off nucleons, nuclei, photons, etc. can be 
legitimately calculated in QCD in the kinematics of fixed x and Q 2 — y oo as 
a consequence of the QCD factorization theorem: tt + T — y 2 jets + T" [62J; 
Jl+N — y V(p, J/*&, p', ■ ■ .)+N' [31], (where the excitation energy of the state 
T' <C Q); j^N — y Meson(7r, K, r], ) + Baryon [56]; Yl ~^ [Few meson system] + 
Baryon [56]; 7* + iV ->■ ^ + iV ^ffl EU E21 ESI EH H OKI S^; and 7* + t -> 
Meson + Meson' [97]. These processes provide new ways to investigate the 
three-dimensional partonic structure of nucleons (transverse distribution of 
partons with a given light-cone fraction) and to compare it to that of A- 
isobars, hyperons, and Ntt. A theoretical analysis of these processes allows 
one also to address such novel questions of short-range parton correlations in 
nucleons as: What is the probability to find a small color singlet cluster in 
the nucleon made of a quark- ant iquark pair, three quarks or even three anti- 
quarks? These processes also probe the minimal light-cone qq components of 
various mesons and few meson systems. In addition, these processes provide 
an effective probe of high energy dynamics of QCD and test whether/at what 
energies the strength of the interaction of small dipoles with nucleons/nuclei 
reaches the maximal strength allowed by unitarity, which leads to breakdown 
of the DGLAP QCD evolution equations. 

An investigations of the same processes off nuclear targets reveals another 
distinctive property of QCD: at fixed x and Q 2 —y 00, nuclear matter is com- 
pletely transparent to the propagation of spatially small colorless clusters of 
quarks and gluons — this regime is usually referred to as color transparency. 
In this limit, the complete transparency for hard diffractive processes unam- 
biguously follows from basic properties of the QCD evolution equations [81] . 
The observation of color transparency of nuclear matter is the striking confir- 
mation that the interaction in QCD is due to the color charge that is screened 
within such clusters. Such a phenomenon would be absent, if high energy 
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processes were dominated by exchanges of usual mesons, as was assumed 
before the emergence of QCD as the theory of hadronic interactions. 

A characteristic feature of these processes is that the final state contains 
a particle (few particles) that has small momentum in the target rest frame. 
Hence to study these processes in the fixed target mode, one needs to design 
a detector that would be able to (1) detect slow particles (including neutrons) 
over a large range of laboratory angles, (2) measure momenta of the leading 
hadrons with high resolution, and (3) operate at high luminosity to reach 
high enough Q 2 — a challenging, though not impossible, task. 

Detection of these reactions in the collider kinematics is somewhat eas- 
ier since the particles that are slow in the target rest frame fly along the 
beam direction. Also, it is much easier to select coherent interactions with 
nuclei. The challenge in the case of these reactions is to reach high enough 
luminosities — so far only channels with vacuum quantum numbers in the 
t-channel were investigated at HERA. 

The prediction and discovery of the quarks-gluon configurations in hadrons 
that weakly interact with a target requires also the presence of configurations 
in the hadron wave function whose interaction with the target is larger than 



average (see the discussion in section 6.4). 



7.2 QCD factorization theorem 

7.2.1 The statement of the theorem 

The starting point for the analysis is the factorization theorem for the process 

ll (<?) + P ^ " Meson" (q + A) + " Baryon" (p - A) (82) 

at large Q 2 , with t and x = Q 2 /(2p ■ q) fixed. It asserts that the amplitude 
has the form of convolution of the three blocks depicted in Fig. [7j 

M = ^ / dz dx 1 f i/p {x 1 ,x 2 ,t,fi)H ij (x 1 /x,Q 2 ,z,fi)(j) j (z,fi) 
i,i Jo J 
+power-suppressed corrections , (83) 

where f/ p is the "generalized parton density" (GPD); X\ — x 2 = x; <j) is the 
light-front wave function of the meson; H ^ is the hard-scattering coefficient 
usefully computable in terms of the powers of a s (Q). The contribution of the 
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Figure 7: The block structure of the DIS exclusive process 72 + p — > 
"meson" + "baryon". 



diagrams, where an extra gluon is exchanged between the hard blocks, is sup- 
pressed by an additional factor of 1/Q 2 - The formal proof [56] is very lengthy 
so we restrict ourselves in the further discussion by the qualitative explana- 
tion only Qualitatively, the factorization in these processes is due to the 
color screening/transparency: the small transverse size of 72 selects small- 
size (point-like) configurations (b ~ 1/Q) in the meson and the interaction 
with such "white" configurations is suppressed by the factor of l/Q 2 - The 
relation of the color screening to factorization is best seen in the Breit frame. 
Before the interaction, 7^ is static, while after the photon is absorbed, the 
quark-gluon system (which would form the meson) moves with a large veloc- 
ity in the direction of the photon keeping the small transverse size, while the 
baryon system rapidly moves in the opposite direction. No soft interactions 
between the left and right movers is possible, provided that the meson sys- 
tem has a small transverse size. The same argument is likely to work for the 
processes where a forward (anti)baryon is produced: 72 +P — >■ forward N + tt, 
1*l + P ~ *" forward A + K + , and 72 + p — > forward p + A^A" [98J, though no 
formal proof has been given so far. 

In the case of the transverse polarization of 7*, the non-perturbative QCD 
contribution is only suppressed by the power of 1/ In Q 2 (similar to the case 
of F2n{x,Q 2 )). It originates from the contribution of highly asymmetric qq 
pairs in the 7^ wave function which have the transverse size similar to that 
of hadrons. 
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7.2.2 Definitions of light-cone distributions and amplitudes: lon- 
gitudinally polarized vector meson 

The wave function of longitudinally polarized vector meson. 
The light-cone wave function of a longitudinally polarized vector meson is 
0JW) = -4= r^e-^"» + (0|V;(y + ,0,O r )7-^(0)|V>,(84) 



where V is a path-ordered exponential of the gluon field along the light-like 
line joining the quark operators entering the matrix element. 

Quark density of the nucleon: 

For a quark of flavor i, its density in the nucleon (qi) reads: 

f l/p (x u x 2 ,t^) = f_%^ e-^ + y-(p'\ T^(0,^,Ot)7 + ^(0) |p)(85) 

Note that in the case of charged mesons, % stands for the flavor indices of the 
initial and final quarks. 

Gluon density of the nucleon: 

For the gluon density in the nucleon, one can give the definition symmetric 
with respect to the X\ — > x 2 transposition: 



fg/ p (X 1 ,X 2 ,t,H) 



dy 1 



( 



~ix2v y 



, 2tt XiX 2 p^ 
x (p'\TG u + (0,y-,O T )VG^(0)\p}. 

(86) 



Note that the factor of l/(x 1 x 2 ) cancels the inverse factor that appears in 
the derivative part of the product of the two gluon field strength tensors 
G u + (0,y~ , Ot)G u+ (0). The normalization condition is 

xf g / p (x, x,t = 0,ti) = fg/ p (x, /i) , (87) 

where f g / p (x,fi) is the usual (diagonal) gluon PDFs. An additional factor of 
x reflects the difference of the symmetric definition of the gluon correlation 
function from that in the diagonal case. 
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The t dependence of the gluon GPDs in the "diagonal" case of X\ = x 2 
is of special importance for the interpretation of various hard pp processes 
(see the discussion in section 7.5). It is described by the normalized two- 
gluon form factor F g (x,t,Q 2 ), where t = — A^_ is the transverse momentum 
transfer to the target. Its Fourier transform describes the transverse spatial 
distribution of gluons with given x: 

F g (x,p\Q 2 ) = ||^ e^) F a (x,t= -AilQ 2 ) , (88) 

where p = \p\ measures the distance from the transverse center of momentum 
of the nucleon. The distribution is normalized such that J cPpF g (x, p\Q 2 ) = 
1. The information on F g (x,t,Q 2 ) that can be extracted from the hard 
exclusive processes like 7 + p — > J/ip + p: 

F g (x,q 2 ,Q 2 ) = 1/(1 -t/m]) 2 (89) 

Here m g is ~ lGeV and slowly decreases with x decrease. 

Modifications necessary for the case of pseudo-scalar meson electropro- 
duction are given in [5*d] . 

The GPDs are different from zero for 

-l<xi<l, -1<x 2 <1. (90) 

There are two physically different regions. In region I, x\ > and x 2 > 0. 
It corresponds to the knockout of a parton with the light-cone fraction x\ of 
the initial target momentum and its absorption in the final state with the 
light-cone fraction x 2 . The Q 2 evolution is described by the DGLAP-type 
evolution equations. In the x — > limit, a simple connection with the diago- 
nal distributions holds. In region II, X\ > and x 2 < 0. This corresponds to 
scattering off a small-size color singlet (qq, gg) emitted by the target. The qq 
case is loosely analogous to scattering off the meson cloud of the target, pro- 
vided the meson is collapsed into a small-size configuration. In this case, the 
Q 2 evolution is similar to the one for the meson wave function and is governed 
by the Brodsky-Lepage-Efremov-Radyushkin evolution equation [55] , 

The imaginary part of the scattering amplitude originates from region I. 
Using a dispersion representation in energy, it is possible to calculate the real 
part of the amplitude for small Xbj and to avoid the consideration of region 
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II. Secondly, at small xy, the space-time evolution of the processes allows for 
a simpler visualization of the interaction process. Also, most of the currently 
available data at large Q 2 were obtained at HERA for the small- x kinematics. 
Thus, as a next step, we summarize the small x-theory and compare it with 
the data. 

7.3 Hard diffractive production of vector mesons 

7.3.1 Space— time evolution of high energy processes 

Vector meson production at small x in the target rest frame can be described 
as a three-stage process [3"3] : 

(i) The longitudinally polarized virtual photon 72 with the four-momentum 
q = (zqo, k t ) breaks up into qq with the lifetime (which follows from the 
energy-time uncertainty principle): 

Ti = lcoh/c = *§— -, ps . (91) 

/Q2 i k ±+ m g m N x 

** "^ z(l-z) 

The coherence length is l co h > 100 fm at HERA. 

(ii) The qq pair then scatters off the target proton, 
(hi) The qq pair then lives for the time 

T* = l f/c = ^ 02) 

z(l-z) 

before the final state vector meson is formed. We note that r/ > r%. 

As a result, the production amplitude A(Yl + P ~ > V + p) can thus be 
written as convolution of the light-cone wave function of the photon, SE^*^^, 
the scattering amplitude of the hadron state, A(nT), and the wave function 
of the vector meson, ipv'- 

A = ¥^ n) (g> A(nT) <g) V(qq -)> V) . (93) 

In the impact parameter space: 

A = Jd 2 b^ rL (b)a(b,s)Mb), (94) 

where b is the transverse separation of q and q. 
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The leading twist expression is 

127r 3 rV_ e+e - M v a 2 s (Q) V 2 v | (l + if ^) xG T (x, Q 2 



rtn L 

ao ~ f *N^VN 



dt 



a EM Q 6 N 2 
t=o 

(95) 

Here, Ty_>. e + e - is the decay width of V — > e + e and 

W ~2 fcted 2 *^^) ( 96 ) 

for Q 2 — > oo. The rapid onset of the leading twist formulas for cr(ee — > 
hadrons) suggests that for p and mesons, $y(z, k t ) and hence 77 are already 
close to the asymptotic value at Q 2 ~ a few GeV 2 . 

Note here that in this expression, the difference between the light-cone 
fractions x\ and X2 was neglected. For large Q 2 , the non-diagonal GPD is 
calculable |100[ 1101] through the diagonal one since the DGLAP evolution 
for GPDs conserves x = X\ — x^, while the light-cone fractions essential at 
the starting point of the evolution grow with an increase of Q 2 . 

In |102j elastic photoproduction of J /if) was evaluated in the leading 
g 2 ln(xo/x) approximation. As we discussed above this approximation ig- 
nores huge NLO effects. 

At extremely small values of x, that are significantly smaller than those 
characteristic for the applicability of the DGLAP approximation, the \b.(xq/x) 
terms not enhanced by \u(Q 2 /Qq) and thus neglected in the DGLAP ap- 
proximation, become important. The restriction on the region of applicabil- 
ity leading log(a;o/a;) approximation follows from the necessity to take into 
account energy-momentum conservation. Indeed, in multi-Regge kinemat- 
ics, the interval in rapidity between adjacent radiations within the ladder is 
Ay 2> 2. This number is comparable with the interval in rapidities achieved 
(to be achieved) in DIS: 

Ay = ln(l/x) + 2 ln(Q/m N ) . (97) 

For the edge of the kinematics achieved at HERA, Ay rj 10. Since four units 
of rapidity are occupied by the two fragmentation regions, two-to-three glu- 
ons are allowed to be radiated in this kinematics. This is insufficient for the 
dominance of multi-Regge kinematics characteristic for LL approximation. 
So far there have been no attempts to describe hard diffractive processes 
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in the resummation approach. Note that formulas obtained in the double 
logarithmic approximation, q s C 1 and a s h\.(xo/x) \n(Q 2 /Ql) ~ 1, should 
coincide for the BFKL and DGLAP approximations. 



7.3.2 Modeling finite-Q 2 effects 



In the convolution integral (Eq. 93 ), \I/f;» (b) with b oc i is convoluted with the 
broad wave function of a light vector mesorQ Hence, the average distances 
contributing to <jl are significantly smaller than those contributing to a?. 
As a result, the effective Q 2 is smaller for vector meson production than for 
&l |103] (Fig. [8]). This effect is taken into account by evaluating cr(qq — N) 
using the dipole model. One also has to include the difference between x\ 
and X2, which is absent when the dipole model is applied to inclusive DIS. 
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Figure 8: The dependence of the average b (left) and the effective Q 2 (right) 
on Q 2 for production of vector mesons |103] . 

A related effect is that, at pre-asymptotic energies, one cannot substitute 
\l/y(6) by \Ev(0). This higher twist correction leads to the suppression of the 
amplitude by the factor: 

T{Q2) = \J d'bdzV^bMqq- N)Mz,b)\ 2 (9g) 

\fd 2 bdz^ rL (z,b)(T(qq-N)0 v (z,O)\ 2 ' 

The HERA data (for a recent summary, see [29]) have confirmed the 
following basic predictions of pQCD 



2 To obtain this expression, one needs first to apply conservation of the electromagnetic 
current to express the bad component of the current through the good one. 
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9 

The rapid increase with energy — \xGn(x, Qlff)\ oc W ' 8 for Q 2 ff ~ 4 
GeV 2 — of p production for Q 2 = 10 — 20 GeV and of J/ip production 
for Q 2 < 10 GeV 2 . (Note that a{W) oc W - 32 for soft physics at t = 
and is even slower for the cross section integrated over t.) For T pro- 
duction, Ql s rj 40 GeV 2 which leads to cxfW) oc VT 17 . This prediction 
maybe tested in the ultraperipheral collisions at the LHC |104j . 

The absolute values of the cross sections of vector meson production 



are well reproduced, provided that the factor T (Eq. 98) is taken into 
account. In the case of T photoproduction, the skewedness effects due 
to large difference between x\ and x 2 as well as the large value of the 
real part of the amplitude are important. Together they increase the 
predicted cross section by a factor of about four |105[ 1106] . 

The decrease of <jl with Q 2 is slower than 1/Q 6 because of the lasC^I 
and T{Q 2 ) factors. 



The ratio cr L /<j T 3> 1 for Q 2 ^> m 



v 



• There is a universal t dependence for large Q 2 originating solely from 
the two-gluon-nucleon form factor. The model, which takes into ac- 
count squeezing of jl with Q 2 , provides a reasonable description of the 
convergence of the t-slopes of light mesons and J/i[) production and 
makes the observation that the slope of J/ip production is practically 
Q 2 independent (Fig. [9]). 

7.3.3 Lessons and open problems 

• Transition from soft to hard regime. 

We can estimate the effective size of a qq dipole as 

B(Q 2 )-B 2g R 2 (dipole) 



B(Q 2 = 0) - B 2g R 2 ' 

where B 2g is the slope of the square of the two-gluon form factor. Based on 
the HERA data [108j . we conclude that 

R 2 (dipole) (Q 2 > 3 GeV 2 )/R 2 p < 1/2 - 1/3 (100) 

for collider energies. Accordingly, it appears that the soft energy dependence 
of the cross section persists over a significant range of the dipole sizes. (This 
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Figure 9: The convergence of the t-slopes, B, of p and J/ip electroproduction 
at high Q 2 . The data are from [1071 1108] ; the curves are the predictions 
of [T03]. 

is consistent with the observed similarity of the energy dependence of p and 
photoproduction.) 

• In the pQCD regime, the t-slope the dipole-nucleon amplitude should 
be a weak function of s = W 2 , B(s) = B(so) + 2a' efi ln(s/so), since the Gribov 
diffusion in the hard regime is small (see the discussion in section 4). Hence, 
a significant contribution to a' comes from the variation of the t dependence 
of the gluon GPD with an decrease in x at Qq. 

• The contribution of soft QCD physics in the overlapping integral be- 
tween the wave functions of the virtual photon and the transversely polar- 
ized vector meson is suppressed by the Sudakov form factor (see, e.g., the 
discussion in |lU9j ). which is absent in the case of the processes initiated by 
longitudinally polarized photons. This is probably relevant for the under- 
standing of similar dependence of ctl,t on x and on t that were observed at 
HERA. 
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7.4 Meson production at intermediate x 

An analysis of the HERA data on vector meson production indicates that, up 
to rather large Q 2 , cross sections of hard exclusive processes are suppressed 
significantly as compared to the leading twist QCD predictions. The origin of 
this is higher twist effects originating from the contribution of the transverse 
sizes that are comparable in the longitudinal photon and the meson wave 
functions. 

At the same time, the overall transverse size of the produced meson is 
quite small (< 0.4 fm) already for Q 2 > 5 GeV 2 . Due to the color trans- 
parency phenomenon, this leads to a strong suppression of the final state 
interaction of the qq pair, which in the end will fragment into the meson and 
the residual baryon system. For W < 20 GeV, this cross section is of the 
order of few mb. Besides, the expansion of the qq system to a normal hadron 
size in the nucleon rest frame takes a significant distance 

hah ~ 2p M /(Am 2 M ) ~ — Lq 2 /(Att4) (101) 

xMn 

where Am 2 M < 1 GeV 2 , and (Am\ I )/2pM is the characteristic light-cone 
energy denominator for a meson M. The condition l co h ^> r^ is satisfied 
for x < 0.2 already for Q 2 > 5 GeV 2 . Hence, it seems likely that the 
precocious factorization into the three blocks (Fig. m) — the overlap integral 
between the photon and the meson wave functions, the hard blob, and the 
skewed distribution — could be valid already at moderately high Q 2 , leading 
to precocious scaling of the spin asymmetries and of the cross section ratios 
as a function of Q 2 . 

The discussion of numerous promising channels such as production of 
charged and neutral vector mesons (p^' , K* + , . . .), pseudoscalar mesons 
7T°, i], rf that are sensitive to the QCD axial anomaly, and A-isobars is beyond 
the scope of this review (for a detailed discussion, see |110j ). 

7.5 Transverse structure of the nucleon at small x from 
GPDs 

At small x many processes are dominated by the two-gluon ladder with no 
contribution from quark GPDs — like production of J ftp and T. In the case 
of GPDs linked to sea quarks, the situation is more complicated. In the case 
of the deeply virtual Compton scattering (DVCS) amplitude in the NLO 

54 



approximation, quark and gluon contributions enter with opposite signs and 
are of a comparable magnitude [lllj . As a result, a relatively small difference 
of the transverse sizes of the sea quarks and gluons, which is expected due to 
pionic cloud effects |112j . is amplified in the t-slope of DVCS. Thus, precision 
measurements of the quark GPDs at small x require an accurate measurement 
of the gluon GPD. 

Higher twist (HT) effects modify the ^-distribution of light mesons up 
to Q 2 ~ 15 GeV 2 . Therefore, it appears that the only practical chance to 
perform a precision measurement of the gluon GPDs is the production of 
onium states. 

The effects of non-diagonality in the gluon GPD appear to be small for the 
J/ijj case, this is due to the large transverse momenta in the wave function, 
which lead to the comparable light-cone fractions of the gluons attached to 
cc: x\jxi ~ 2, x\ — %2 = x = (Q 2 + m 2 ,,)/W 2 . In the T case, the effect of 
non-diagonality can be taken into account via the DGLAP evolution. As a 
result, 

ex F 2 (x,t)exp(ABt) , (102) 

where F g (x, t) is the two-gluon form factor of the nucleon; the second factor 
takes into account a small but finite correction due to the finite size of J /if} 
that was estimated in |lU3j to be AB rj 0.3 GeV -2 . 

The ^-dependence of the measured differential cross sections of exclusive 
processes for \t\ < 1 GeV 2 is commonly described either by an exponential or 
by a dipole form inspired by analogy with the nucleon elastic form factors. 
The data are not precise enough to distinguish between the two forms since 
they mostly differ at small t, where the resolution of the measurements is 
moderate, and at large —t > 0.8 GeV 2 , where the measurements rather 
strongly depend on the procedure of subtraction of the inelastic background 
contribution. 

The data can be fitted in the following form [113J: 

B g (x) = B g0 + 2a' g lm>o/x), (103) 

where 
a* = 0.001, B g0 = 4.1 (t° j) GeV~ 2 , a' g = 0.140 (i°;°g) GeV" 2 . (104) 
Fits of similar quality are produced with the dipole form: 

F g (x,t\Q 2 ) = (1-t/m 2 )" 2 , B g = 3.2/m 2 g (form 2 - 1 GeV 2 ). (105) 
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The spatial distributions of gluons in the transverse plane corresponding to 
the two fits are: 

f (2nB g )- 1 exp[-p 2 /(25,)] , 
F g (x,p\Q 2 ) = \ (106) 

{ [m 2 g /(2n)}(m g p/2)K 1 (m g p), 

These transverse distributions are similar for the average p, leading, for ex- 
ample, to the nearly identical distributions over the impact parameter for 
production of dijets in pp collisions at the LHC |113j . At the same time, 
the dipole fit gives a significantly larger F g (x, p\Q 2 ) for small b. As a result, 
analyses of the proximity to the black disc regime for the interaction of a 
small dipole with the nucleon for small b are sensitive to the choice of the 
model for F g (x,p\Q 2 ). A related effect is a factor of 1.6 difference of lb 2 ') 
in the two fits. The conclusion |114j that gluons are localized in a smaller 
transverse area than that given by the e.m. form factor is based on the use 
of the same shape for the e.m. and two-gluon form factors. 

The current knowledge of F g (x,p\Q 2 ) allows us to study the impact pa- 
rameter dependence of dijet production at collider energies. One finds that 
the median impact parameters for the inelastic pp collisions with jet pro- 
duction are a factor of two smaller than those for the minimal bias inelastic 
events, and weakly depend on the rapidity and pt of the jets. This may 
explain the regularities in the multiplicity of the underlying events on p t of 
the trigger |113| . It also provides an important constraint for the models 
of inelastic pp collisions at the LHC and, in particular, of the dynamics of 
multiparton interactions. 

7.6 Break-up processes with gaps: from small to large t 

7.6.1 Probing fluctuations of the gluon field 

In the high energy 7* + p — > M + rapidity gap + X process, at £ ~ the 
two-gluon ladder couples only to one parton in the target in the leading 
twist approximation. If the strength of the coupling to all configurations 
containing partons with a given x were the same, it would be impossible to 
produce an inelastic final system X. As a result, similarly to the case of 
inelastic diffraction of hadrons off hadrons, the discussed process measures 
the variance of the gluon field. It is given by the ratio of the diffraction 
dissociation and elastic cross sections for vector meson production at t — 
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Our estimates of the strength of the fluctuations due to the fluctuations of 
the overall size of the nucleon, as seen in soft inelastic diffraction, find that 



Uhard 



0.15 at small x, which is consistent with the current data. 




Figure 10: A typical QCD diagram for the rapidity gap process (108). 



7.6.2 Onset of new regime at large t 

In order to determine the spacial distribution of gluons at b ~ 0, one needs 
to study exclusive processes at large t since in this case the integral over 
—t in the Fourier transform converges very slowly. For example, in the case 
of the dipole fit, half of the contribution to F g (x,b = 0) originates from 
—t > m? and one quarter, from —t > 3m; ?. This implies that one needs a 
detector which will be able to separate exclusive processes from diffractive 
dissociation in a very broad range of t, where the dissociation dominates by 
far. 

A further complication is that the QCD factorization theorem [56] was 
derived in the limit of t = const, x = const, and Q 2 — >■ oo. Recent studies [33] 
found that the pattern of the QCD evolution changes in the kinematical 
domain when —t becomes comparable to the intrinsic hardness scale of the 
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process {Q 2 S ~ 3 GeV 2 for J /if) photoproduction). The DGLAP equations 
require modifications in this case. Let us, for example, consider the process 

7* + p -> J/ip + rapidity gap + M x (108) 

for large — t, which is still smaller than Q 2 + M 2 T/ ,. The typical leading QCD 



L j/i>- 



diagrams (Fig. 10) correspond to the attachment of the two-gluon ladder to a 
parton with xj = —t/(—t + M\ — m 2 N ). The cross section has the factorized 
form similar to [116J: 

UO7+P— >V+X QC"7+quark— ^-V+quark 



dtdxj dt 



Q1 

—g P (x,j, t) + ^2(q l p (xj, t) + c? p (xj, £)) 

(109) 

At fixed Xj, the energy dependence is determined by the evolution of the 

dipole-parton elastic scattering amplitude with x/xj. In difference from the 

small t limit, the DGLAP evolution is strongly suppressed and completely 

disappears for —t close to the intrinsic scale. The HERA data J117] on 



the energy dependence of the process (108) are consistent with the behavior 
expected in QCD. 

Consequently, the effective ajp(t) in this limit stays close to unity until 
very small x/xj (not available at HERA) where the Pomeron-type behavior 
may reveal itself. 

The same mechanism may be responsible for part of the drop of ajp(t) 
with an increase in —t observed for elastic J/ip production. This may indi- 
cate that precision measurements of the t dependence of the two-gluon form 
factor at high energies would require using of electroproduction rather than 
of photoproduction. The discussed phenomenon may also be relevant for the 
explanation of the pattern observed in photoproduction of p mesons, where 
ajp(t) appears to flatten out around ajp(t) = 1 for large \t\ (Fig. [3 

7.6.3 Probing minimal quark component in the pion 

The QCD factorization theorem discussed in the previous subsections allows 
us to calculate another group of hard processes where the selection of the final 
state dictates squeezing of the initial state. This is a particular case of the 
pre-selection phenomenon familiar from non-relativistic quantum mechanics. 
The most straightforward process is 

7r + T — > two jets + rapidity gap + T . (HO) 
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This process is in a sense a mirror image of vector meson production in DIS. 
The pion in the initial state collapses into a small-size configuration due to 
the hard interaction; this qq pair interacts coherently with the target and 



transforms into two jets [62] (Fig. 11) 




Z, K t 



I -Z, -Kt 



A(N) ^B^— A(N) 



Figure 11: A typical two-gluon ladder exchange diagram contributing to the 
process of pion coherent diffractive dissociation. 

In the limit of large transverse momenta of the jets, one can justify ap- 
plicability of the QCD factorization theorem and obtain [118J: 

da(ir + A^2jet + A)(q t = 0) _ (1 + r? 2 ) \ « s tt 2 2 ] 2 

dFdz~d^ t " 16^(2^ [ A {XAZ > Kt)) ~^ X ^ A ^ x ^ Q ) 

(HI) 
where x^{ z -, K t) = 4ttCf * \3fnz(l- — z)', A is the Laplacian in the n t space; 
z is the fraction of the light-cone plus-momentum carried by the quark in the 
final state; XiGa{xi,X2, k 2 ) is the generalized gluon density of the nucleus, 
where X\ and X2 are the fractions of the target momentum carried by ex- 
changed gluons 1 and 2 respectively, X\ — x 2 = M^- et /s and x 2 < X\ and 
the integral over x 2 is not written explicitly; rj = ReF/ImF, where F as the 
dipole-nucleon scattering amplitude). Note that the resulting n~j~ 8 depen- 
dence is a consequence of the comparatively well-understood wave function 
of the pion. This wave function determines the asymptotic behavior of the 
pion electromagnetic form factor. 

For x > 0.03, Ga(x, Q 2 ) = AGn(x,Q 2 ). Hence, for this kinematics there 
are no absorptive effects — the amplitude at t = should be oc A. This predic- 
tion has been confirmed at FNAL with p„.=500 GeV/c: the experiment |119j 



59 



observed a strong coherent peak for dijet production from carbon and plat- 
inum targets and measured the A dependence for this interval of A that was 
found to be A 1,54 . The ratio of the cross sections for the two targets is a 
factor of seven larger than that for soft coherent diffraction. Furthermore, 
the observed dependence of the cross section on the pion momentum fraction 
and the jet transverse momentum is well-consistent with the perturbative 
QCD prediction of [62 1 1118] for kt > 1.5 GeV/c. The relatively early onset 
of scaling for this process — as compared to diffract ive electroproduction of 
vector mesons discussed above — maybe due to the presence of the plane wave 
in the convolution formula for the dijet production cross section. This should 
be compared to the case of the virtual photon wave function that restricts 
the phase space much stronger. 

We would like to note that the derivation of the QCD factorization the- 
orem for this process |118j heavily used the fact that the trigger on two 
high-K t jets along with the Ward identities enforce the qq pair and, there- 
fore, the color to be concentrated in the interaction volume oc 1/ftf. The 
use of the asymptotic freedom and the restriction by the leading twist con- 
tribution guarantees the dominance of the gg-component in the pion wave 
function at sufficiently large k. The generalized gluon distribution of the 
target depends on variables Xi 3> x 2 and X\ — x 2 = M|j et /s. The kinemat- 
ics follows from energy-momentum conservation. The mass of the two jets, 
M 2 j et , is significantly larger than the pion mass; that allowed one to justify 
the applicability of pQCD and the Ward identities. The calculation also ex- 
plores the asymptotic solution of the QCD evolution equation for the pion 
wave function. 

It was assumed in [120j that for the calculation, one can substitute the 
pion by a system of the non-interacting quark and antiquark with the mass 
equal that of the two-jet system M^et 3> m n . This assumption violates the 
conditions of the applicability of the QCD factorization theorem discussed 
above and, therefore, produces a different amplitude. There is no requirement 
for the qq pair (color) to be within a small volume. Therefore, there is 
no way to justify the applicability of pQCD and to neglect other quark- 
gluon configurations. The calculation of leading Feynman diagrams due to 
the two-gluon exchange produces the factor l/(M|= et — M|), while in the 
kinematics of [120j, Mi = M 2 - ]e t- A prescription is needed to remove this 
artificial singularity. Violation of the Ward identities requires an additional 
prescription of how to derive the condition that the leading amplitude for 
qq scattering off the two-gluon ladder exchange is oc rf, where r t is the 
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momentum of the exchanged gluon. In this kinematics, X\ = x<i ~ 0. Besides, 
in the framework of the approximations discussed above, the authors |120j 
found that the scattering amplitude contains singularities that contradict the 
Landau rules for the amplitudes of physical processes. 

For a proton projectile, the related process would be proton diffraction 
into three jets: 

p + A — > three jets + rapidity gap + A . (H2) 

So far, the observation of this process at the LHC looks very difficult — current 
detectors do not have acceptance for realistic p t of jets of this process. 

7.7 Nuclear effects in hard exclusive reactions 

The use of nuclear targets provides complementary probes of the QCD dy- 
namics of diffractive processes. QCD predicts that a sufficiently energetic, 
spatially small color-neutral wave packet of quarks and gluons should tra- 
verse hadron medium without absorption. This prediction is equivalent to 
the QCD factorization theorem of [SB]. The complete transparency of nu- 
clear matter [color transparency (CT)] has been unambiguously observed at 
several experiments at high energies (see discussion below). Such processes 
are becoming a promising tool for the detailed investigation of the quark and 
gluon structure of nuclei. 

We give here just a few examples. The factorization theorem [56J predicts 
for coherent electroproduction of vector mesons 



da 
dt 



(YA^VA)\ t __ 



da 
dt 



(j*N^VN)\ t __ 



Fk(x,Qy :l 

Fk(x,Q) 



G a( X iQ) _ A2a g (x,Q) t U o\ 

G%(x,Q)~ A ' (U3) 



with a g (x, Q) ~ 1 for x > 0.02. Thus final-state interaction in this process 
is a higher twist effect. 

• Almost complete transparency has been observed at FNAL |121j in both 
the coherent and incoherent processes: 'jA — > J/\l/ + A, 7 A — > J/\l/ + A*. 

• Complete transparency of nuclear matter, which follows from the QCD 
factorization theorem for the sufficiently large transverse momenta of jets 
kt > 2 GeV, has been observed in the processes 7r + A — > 2 jets + A' [119J. 

• For quasi-elastic scattering off nuclei and for large enough |t| > 0.1 
GeV 2 , one can use closure over the processes of nuclear disintegration A'. 
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For example, for production of neutral mesons, one obtains 



£ 



da{Y L + A^M + A') _ da{Y L + P -> M + p') 

Zj 



dt dt 



+ N M?L + n^M + n>) ^ (n4) 

( / f 



• The color transparency phenomenon leads to a strong suppression of 
the cross section of coherent scattering off the lightest nuclei for the t range, 
where double scattering dominates |122j . 

At lower energies, where the Lorentz factor is not large enough to guar- 
antee the sufficiently large lifetime of a spatially small wave packet of quarks 
and gluons traversing a nuclear target, CT is masked to large extent by the 
following quantum mechanical phenomenon: the small-size wave packet is not 
an eigenstate of the QCD Hamiltonian and, hence, it rapidly expands [123J. 

• The onset of complete transparency has been observed |124] in the in- 
coherent process 7*+A->p + Xati/~ 200 GeV in the kinematics 
where hadron production in the nuclear fragmentation region was al- 
lowed. A smaller effect was observed at v ~ 15 GeV under similar 
conditions [125J. 

• A gradual increase of transparency in exclusive processes with an in- 
crease of Q 2 was observed in high precision experiments at JLab at 
v ~ a few GeV that measured production of 7r + |126] and p° [127] . 
The observed A and Q 2 dependencies are consistent with the familiar 
quantum mechanical effect of the expansion of a small qq wave packet 
with distance 



In collider kinematics, it is much easier to check whether the nucleus 
remained intact and to also select certain exclusive break-up channels. Hence, 
it may be easier to study color transparency effects for these processes at 
colliders. 

Overall, the use of nuclei would add substantially to the program of stud- 
ies on hard exclusive processes for the x > 0.05 range. Since nuclear parton 
densities for a wide range of x > 0.05 are oc A up to a small correction due 
to the EMC effect, one would not have to deal simultaneously with the effect 
of leading twist shadowing. 
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7.8 Summary 

Hard exclusive meson (few meson) production is calculable in QCD in the 
same sense as the leading twist DIS processes. Hard exclusive processes 
provide unique ways to study minimal Fock state components of mesons and 
structure functions of mesons, to compare skewed parton distributions in a 
multitude of baryons, and to investigate the onset of color transparency. 

Colliders have a number of advantages for observing many of these pro- 
cesses. 

To have a successful program of studies on the hadron and nuclear struc- 
ture in the next decade, one needs a triad of electron accelerator facilities 
optimized to study the following three x ranges: (a) the x > 0.3 range rele- 
vant for the study of short-range correlations in nucleons and nuclei, (b) the 
0.05 < x < 0.3 range relevant for the study of multiparton correlations in 
nucleons and the origin of the nuclear forces, (c) the x < 0.05 range rele- 
vant for the nuclear shadowing phenomena, high parton density physics, etc. 
This corresponds well to the W < 8 GeV range discussed for the upgrade of 
the Jefferson Lab facility and for the COMPASS experiment with the recoil 
detector as well as to the W ~ 30 — 150 GeV range discussed for the electron- 
nucleon/nucleus colliders at JLab and RHIC and for the LHeC collider at 
CERN. 



8 A new regime of high energy QCD 

8.1 Introduction 

At large energies, hard QCD interactions become strong at central impact 
parameters in spite of the small running coupling constant and a larger scale 
of momenta. The necessity to take into account an entire series in (1/Q 2 ) n 
(which resembles the situation in the second-order phase transitions) shows 
that this new phase of QCD has a different continuous symmetry (conformal 
symmetry), thereby distinguishing it from the perturbative phase of QCD. 

With a further increase of energy, predictions based on the QCD factor- 
ization theorems start to contradict probability conservation. For the inter- 
action of a spatially small dipole with a target, this can be be formulated 
as: 

a e i(dipole + T) < -<j tot (dipole + T) . (H5) 
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The upper limit on the (T e \/cr tot ratio is reached in the regime of complete 
absorption for all essential impact parameters. The restriction due to proba- 
bility conservation is stronger for the scattering at central impact parameters 



(see the discussion in section 8.2). 



In the leading twist approximation, the inequality (115) is violated at 
sufficiently small x since the QCD factorization theorem for the interaction 
of a small-size colorless dipole leads to cr c i oc (xGt(x,Q 2 )/Q 2 ) 2 /B, which 
becomes larger than er to t ^ xGt(x,Q 2 )/Q 2 for sufficiently small x (B is the 
t-slope of the differential cross section of the dipole-target elastic scattering). 
Thus, in the limit of fixed and large Q 2 , x — > and fixed impact parame- 
ters, the decomposition over high powers of {1/Q 2 ) n (over twists) becomes 
ineffective. Indeed, the dependence of the higher twist term T n+ 2 on x and 
Q 2 can be easily evaluated at large energies: T n+2 oc (1 / Q 2 ) n (xq / xY 1+n ^ x ^ 1 \ 
Therefore the ranking over twists disappears at sufficiently small i. Asa re- 
sult, the concept of a spatially small dipole becomes ineffective as well since 
the contribution of various higher twist effects (for example, the splitting of 
a small dipole into two small dipoles each interacting with the target) is not 
suppressed at these energies. In other words, the effective number of dipoles 
continuously increases with energy. 

The regime of complete absorption at fixed impact parameters and con- 
servation of probability does not preclude a rapid increase in cross sections 
of hard processes in pp collisions at central impact parameters as well as of 
the 7*iV cross section with a decrease in x at fixed Q 2 . At collider energies, 
the hard contribution to the nucleon structure functions increases with en- 
ergy as 0"(7*iV) oc ln 2 (xo/x), which is faster than the x~ ' 2 behavior observed 
at HERA at Q 2 ~ a few GeV 2 . At very high energies, an even faster in- 
crease is expected, <r(7*iV) ex ln 3 (x /a;) \V2S\ 1129] . The additional ln(x /x) 
is a consequence of the singular behavior of the light-cone wave function of 
the virtual photon in the coordinate space, which in the momentum space 
corresponds to the contribution of the quarks with momenta ^> Q in the 
box diagram. Numerical studies on the energy dependence of F 2 n(x,Q 2 ) 
that took into account the taming of partial waves at small impact parame- 
ters, were performed in a number of papers using the dipole model (see, for 
example, [130J). 

Note also that QCD predicts taming of structure functions of nuclei which 
competes with the nuclear shadowing phenomenon. However, it is hard to 
observe the violation of the DGLAP approximation by an analysis of the 
experimental data in the current energy range since the evolution equation 
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is linear and has the flexible initial condition. 

Deep inelastic interactions studied so far at HERA are far from the strong 
absorption regime for the QCD interactions with a small coupling constant. 
The only possible exception is the the gluon dipole-nucleon interaction at 
central impact parameters. However, the situation may change in the case of 
scattering off heavy nuclei and also in the case of central pp/pA/AA collisions 
at the LHC. 

8.2 Conflict of pQCD with probability conservation. 

Asymptotic freedom in pQCD does not guarantee probability conservation 
since pQCD amplitudes are rapidly increasing with energy. The conflict with 
probability conservation reveals itself in the scattering of a color singlet wave 
packet of quarks and gluons of the small diameter doff a hadron target, see 
Eq. 115| The cross section of elastic scattering of a dipole off the nucleon 



target can be parametrized as: 

4- /\H rr, .(rlinnlp 4- N\ 2 

Fl(t), (116) 



da e \(dipole + N — > dipole + N) a to t(dipole + N)' 2 
It ~ 16tt 2ffV 



where F2 9 (t) is the two-gluon form factor of the nucleon extracted in |114j 
from the HERA data on hard diffractive electroproduction of J/ip- At t = 0, 



Eq. 116 is the optical theorem. In the expression above we neglected a small 



contribution of the real part of the elastic dipole scattering amplitude. Using 



Eqs. |116| and |115[ we obtain the following inequality: 
cr tot (dipole + N) 



16vr 



JdtF* g (t)< 1/2. (117) 



In the case of scattering at the zero impact parameter, the inequality is a 
factor of ~ 2 stronger: 



o- to t(dipole + N) 



JdtF 2g {t)<\. (118) 



87T 

One can also study the unitarity condition as a function of an impact pa- 



rameter using Eq. 57 |131[ 1132] . 

The cross section for a spatially small color singlet dipole scattering off 
the nucleon target follows directly from the pQCD factorization theorem 
in the leading a s hi{Q 2 /Q"fy approximation [IJU E21 E3]- It can also be de- 
rived from the Born term obtained in the leading a s ln(xo/x) approximation 
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in [61] . Using the machinery of pQCD calculations, it should not be too dif- 
ficult to calculate a series of LO, NLO, and NNLO approximations including 
corrections to the cross section and to the dipole wave function itself. 

Since the expressions in the leading log approximation are too bulky, for 
illustration purposes, we use the fit xG oc (x /x) x ^ Q '. Due to an increase 
in X(Q 2 ) with Q 2 , the drop of the amplitude due to a decrease in the dipole 
size is compensated to some extent by a faster increase in the gluon density. 
Thus, x cr at which unitarity breaks down, increases with an increase in Q 2 
rather slowly. 

A similar conflict with probability conservation exists for the perturba- 
tive " Pomeron" . The energy at which this conflict becomes acute obviously 
depends on the impact parameter, see |130] and references therein. Note also 
that the conflict with unitary exists not only for scattering of small dipoles 
off hadrons, but also in the case of scattering of two small objects [133J. 

8.3 Gluon dipole— proton scattering at the upper edge 
of HERA kinematics 

To demonstrate that the new QCD regime for a gluon cloud is not far from 
the kinematical domain probed at the LHC, we consider here scattering 
of a colorless small gluon dipole off the nucleon (for a detailed discussion, 
see |134j ). One can evaluate Y gg (x, b) based on the information on the total 
cross section of the qq dipole-nucleon interaction extracted from DIS inclu- 
sive data (Eq. 65), the t dependence of the two-gluon form factor (section [7]) , 



and the relation 



q 

-unel/„, i,\ -pincl/ 



r r (*> b ) = -f7i (*> h ) ■ ( 119 ) 

The knowledge of T gg (x,b) allows us to calculate the total and elastic cross 
sections of the colorless " gluon-gluon" dipole-nucleon scattering. (By defini- 
tion, the imaginary part of the partial wave for the dipole-nucleon scattering 
is T gg (x, b).) The ratio of the elastic and total cross sections for the scattering 
of a color-singlet dipole off the nucleon is: 

, ^ 2N [<fb\r q Jb)\ 2 , N 

Numerical estimates |134] indicate that for d = 0.1 fm (Q 2 ~ 40 GeV 2 ) 
corresponding to T photoproduction, the value of R g (x, Q 2 ) is small (~ 0.14) 
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for x ~ 10~ 4 . However, R g (x,Q 2 ) reaches the value of ~ 0.4, when x goes 
down to 10~ 7 indicating that for such x, most of the cross section maybe 
due to the interaction in the BDR. At the same time, for Q 2 ~ 4 GeV 2 one 
finds that the interaction is nearly black over a rather large range of impact 
parameters already for x ~ 10~ 4 , leading to R g ~ 0.3 — 0.4. At the same 
time, in the gg-nucleon case, significantly smaller x are necessary to reach 
the BDR. So the interaction in this channel is rather far from the regime of 
complete absorption at x achieved at HERA. 

The predictions discussed above can be compared with the experimental 
data on diffraction in DIS obtained at HERA. The sum of the cross sections 
of color-singlet gluon dipole elastic scattering and of inelastic diffraction scat- 
tering is calculable in terms of diffractive gluon PDFs (when an external hard 
probe couples directly to the gluon) integrated over the momentum of the 
diffracted proton with the cut xp < 0.03. The total cross section of the color 
singlet dipole-proton scattering is calculable in terms of the gluon density of 
the nucleon. Thus, to compare R g with the data, one should measure the 
ratio of the diffractive and inclusive cross sections induced by a hard probe 
coupled to gluons: 

R / Q2) = L 0,03 dx P Jo"" dtf? (4) (x/x P , Q 2 , x F , t) 

where f g is the diffractive gluon PDF. It can be evaluated using the 
most recent analyses of hard diffraction at HERA. One finds that R g (x ~ 
1(T 4 , Q 2 = 4 GeV 2 ) ~ 0.3 and R q (x ~ 10" 4 , Q 2 = 4 GeV 2 ) ~ 0.2 which 
confirms proximity of the gluon interaction to the BDR at small impact pa- 
rameters (for the recent update, see Fig. 69 in [135] ). 

A word of caution is necessary here. These calculations are based on 
the leading log approximation where such higher order Fock states as qqg 
are neglected. Hence, these calculation should only be considered as semi- 
quantitative ones. 

8.4 The regime of complete absorption 

Bjorken scaling, i.e., the dependence of cross sections of hard processes only 
on the hard scale, completely disappears in the black-disc regime only for 
a specific range of Q 2 that depends on x. Basic formulas of the black-disc 
regime for the total cross section of photo(electro)production of hadrons off 
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a heavy nucleus at high energies were derived by V. Gribov [136J by analyz- 
ing the contribution of diffraction to the structure functions of heavy nuclei 
and an increase in absorption with an increase in the atomic number. The 
theoretical observation that enabled the calculation was that in the BDR, 
non-diagonal transitions between diffractively produced hadronic states are 
absent at zero angles. 

Interesting features of this regime include the dominance of diffraction 
into dijets which constitutes 50% of the total cross section, a gross change 
of the Q 2 dependence of exclusive meson production (from 1/Q 6 to 1/Q 2 for 
the case of the longitudinal photon), and gross suppression of leading hadron 
production in the current fragmentation region (suppression of the effective 
fractional energy losses) |128j . 

A big challenge is to establish the mechanism of the onset of the BDR 
and the basic features of the BDR. A popular approach is to tame the rapid 
increase in the perturbative LO "Pomeron" with energy by taking into ac- 
count the simplest non-linear effects [137J. This approach has even more se- 
vere problems as compared to the perturbative " Pomeron" . Indeed, energy- 
momentum conservation restricts the number of allowed branchings of the 
"Pomerons" due to the triple "Pomeron" and the number of mult i-" Pomeron" 
exchanges. The restriction follows from the fact that an exchange of the per- 
turbative "Pomeron" is dominated by inelastic processes which carry energy 
and momentum. This restriction is especially severe for the multi-Regge kine- 
matics explored in |137] . Thus, the contribution of fan diagrams relevant in 
this approach for the screening of the perturbative " Pomeron" requires colli- 
sion energies that cannot be treated in a realistic way at the collider energies. 

The appearance of new scales in addition to Q 2 follows from the fact that 
the decomposition over twists becomes meaningless at sufficiently small x 
since the decomposition parameter starts to exceed the radius of convergence. 
As a result, an analytic continuation becomes necessary. 

The new scale depends on the incident energy as a result of probability 
conservation. This property of the new QCD regime distinguishes it from 
the pQCD regime. It may indicate that in this regime, violation of the two- 
dimensional translation symmetry, or even conformal invariance, takes place. 
This option has been suggested for soft processes in [1381 1139] and for hard 
QCD in |129j . In this case, the quasi-Goldstone mode is the motion along 
the impact parameter in the region of complete absorption. Moreover it is 
unclear whether pQCD is applicable even qualitatively for the description of 
the new QCD regime since in the kinematics where the interaction becomes 
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strong, one should take into account the Gribov copies found in [140] . 

9 Conclusions 

The theory of the S'-matrix allowed one to understand many regularities of 
diffractive processes in the regime that is now referred to as the soft QCD 
regime. Further progress became possible with an advent of QCD and the 
focus on the processes involving a hard scale. 

Investigations of diffractive phenomena led to the ideas, concepts, and 
technologies of the calculations that form the foundation of modern particle 
physics and are the basis for the future developments aiming at bridging 
the gap between soft and hard phenomena. The exploration of the color 
transparency phenomenon supplies a new method of investigations of hadrons 
and nuclei. Investigation of the new QCD regime of strong interaction will 
allow one to find new phenomena and to develop new methods of treating 
phase transitions in the relativistic kinematics. The long-standing challenges 
include confinement of quarks and gluons and spontaneously broken chiral 
symmetry. 

Also, QCD is the only quantum field theory leading to non-linear phe- 
nomena that can be probed in laboratory. Hence, it may provide new tools 
for developing a theory beyond the Standard Model of strong and electroweak 
interactions. 
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